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Abstract. In a previous paper by the first two authors, a tube formula for fractal sprays 
was obtained which also applies to a certain class of self-similar fractals. The proof of 
this formula uses distributional techniques and requires fairly strong conditions on the 
geometry of the tiling (specifically, the inner tube foriTiula for each generator of the fractal 
spray is required to be polynomial). Now we extend and strengthen the tube formula by 
removing the conditions on the geometry of the generators, and also by giving a proof 
which holds pointwise, rather than distributionally. Hence, our results for fractal sprays 
extend to higher dimensions the pointwise tube formula for (1 -dimensional) fractal strings 
obtained earlier by Lapidus and van Frankenhuijsen. 

Our pointwise tube formulas are expressed as a sum of the residues of the "tubular 
zeta function" of the fractal spray in R"*. This sum ranges over the complex diiTiensions 
of the spray, that is, over the poles of the geometric zeta function of the underlying fractal 
string and the integers Q,i,. . .,d. The resulting "fractal tube formulas" are applied to the 
important special case of self-similar tilings, but are also illustrated in other geometrically 
natural situations. Our tube fomiulas may also be seen as fractal analogues of the classical 
Steiner formula. 
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1. Introduction 

Our main results are tube formulas for fractal sprays in R'', the higher-dimensional 
analogues of (geometric) fractal strings in R. In [LvF2], a fractal string is defined to be a 
bounded open subset of the real line R; see also [Lapl, Lap2, Lap3, LaMa, LaPol, LaPo2, 
LvF3, LvFl, LaPel]. Here, we emphasize the interpretation of a fractal string as a sequence 
of positive numbers, rather than as a collection of open intervals in the geometric sense. 

Definition 1.1 (Fractal string). A fractal string £, - is a nonincreasing sequence of 

positive real numbers Ij > satisfying limy^co (j - 0. 

In particular, we do not assume that the sum of the lengths of a fractal string is finite. 
Hence, X. might not have a geometric realization as a bounded open subset of R, in the 
sense of [LvF2, Def. 1.2]. 

Definition 1.2 (Fractal spray). A fractal spray T defined on a bounded open set U Q W' 
via the fractal string X - K^l^li is ^ collection of disjoint bounded open sets {U^}J^^ in 
R'' such that, for each set with j > 1, there exists a similarity transformation of R'' 
with scaling ratio fj and satisfying f/' = '¥j(U). The spray T is said to be scaled by the 
fractal string X, and the connected components of the set U are called the generators of 
the fractal spray. The generators are denoted by Gg, where q ranges over some finite or 
countable index set. When there is only one generator, we denote it by G instead of Gi . 

Hence, a fractal spray on the generator G is just a collection of disjoint scaled copies of 
G such that the scaling ratios form a fractal string (in the sense of Definition 1 . 1 ), just as in 
[LaPo2] and [LvF2, Section 1.4]. Note that since U is bounded and open, each generator 
is a bounded open and connected subset of R'', and hence there can be at most countably 
many generators. We always assume in the sequel that T has finitely many generators 
{G^l^^j, which allows us to study only the case of a single generator G (see the explanation 
at the start of Section 4, and the discussion just following (5.7)). 

Note that fractal strings in the geometric sense may be viewed as fractal sprays in R 
generated by a bounded open interval G; indeed, they are disjoint unions of a sequence 
of bounded open intervals. Therefore, geometric fractal strings are included in the setting 
of fractal sprays. An important subclass of fractal sprays is formed by self-similar tilings, 
which appear naturally in connection with self-similar sets and are higher-dimensional 
generalizations of the (geometric) self-similar strings studied in [LvF3, LvFl, LvF2]; see 
Section 5. 

In the classical literature, the s-parallel set (or s-neighborhood) of a bounded set A c R'' 
is the set of points within (Euclidean) distance e of A (see (1.5)), and a tube formula for A 
is an explicit expression for the volume of the e-parallel set of A, viewed as a function of 
e; see Section 1.1. In this paper, we make use of the following "inner" analogues of these 
notions: 

For e > 0, the inner e-parallel set (or inner s-neighborhood) of a bounded open set 
A c R'' is the set 

A_s:={;ceA -. dist(x,A'^)<e), (1.1) 

and an (inner) tube formula for A is an expression giving the volume V{A,e) :- /{^(A-e) 
(i.e., the (^-dimensional Lebesgue measure) of the set A_£ as a function of e € (0, oo). 

Similarly, by a tube formula for a fractal spray T - {G^'J^y, we will simply understand 
an expression V{T, s) for the volume of the inner e-parallel set r_£ of the union set T := 
U^i of the components G^ as a function of e; that is. 
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V(T,e) - AAT^,) = ^^^((GO-.) = J] ^(^''^)- 



(1.2) 



Our main results in this paper are tube formulas for fractal sprays in R'', which are given 
in Theorem 4. 1 and Corollary 4.2, and later specialized to the class of self-similar tilings in 
Theorems 5.7 and 5.12, along with their respective corollaries, the 'fractal tube formulas' 
obtained in Corollaries 5.9-5.10 and Corollary 5.13. 

These tube formulas express the volume V(T, s) of the inner e-parallel sets r_£ of the 
given fractal spray T in R"' with d > 1, as a (typically infinite) sum over the set Dt £ C 
of complex dimensions of T. These complex dimensions are defined as the poles of the 
tubular zeta function = ^r(^' associated with the spray T (see Definition 3.5 and 
Proposition 3.6), and each summand is equal to the residue' of at the corresponding 
complex dimension. Roughly speaking, we show that for all sufficiently small e > 0, 



where is suitably defined in terms of the scaling zeta function ^ of the underlying fractal 
string X. and the geometry of the generator G of the spray. Here, ^(i) is the meromorphic 
continuation of the Dirichlet series 2" j ^J, initially defined for Re(i) sufficiently large; 
see Definition 3.1. Moreover, the set Dt = 2)£ U {0,1, . . . ,d} of complex dimensions of 
T consists of the integer dimensions 0,1, ... ,d and the scaling dimensions, which are the 
complex dimensions of the associated fractal string X. (defined as the poles of the scaling 
zeta function {£); see Definitions 3.7 and 3.1. 

Depending on the assumptions, our pointwise tube formulas are either exact (as in ( 1 .3) 
just above) or else contain an error term, which can be estimated explicitly as s tends to 
zero. In the latter case, the aforementioned sum of residues ranges only over the 'visible' 
complex dimensions of T; i.e., those complex dimensions lying in a window W, the region 
to the right of a (suitably chosen) vertical contour S called the screen; see Section 3 for the 
precise definitions. 

The fractal tube formulas obtained in this paper extend previous results in two ways. 
First, we extend the scope of the tube formulas of [LaPe3] to fractal sprays whose genera- 
tors may be arbitrary bounded open sets. Second, we give a pointwise version of the tube 
formula obtained in a distributional sense in [LaPe3]. This generalizes and clarifies the 
results previously obtained in [LaPe3, LaPe2, PW, Pea, LvF2, LaPel]. 

Furthermore, formula (1.3) (given precisely in Theorem 4.1) directly extends [LvF2, 
Thm. 8.7] (the pointwise tube formula for fractal strings) to higher dimensions and implies 
as a corollary that all self-similar tilings have a fractal tube formula; see Section 4.1 and 
Section 5 for more details. The tube formulas with error term also allow us to obtain 
information concerning the Minkowski measurability (and Minkowski content, when it 
exists) of fractal sprays and self-similar tilings under certain conditions; this is taken up in 
[LPWl]. 

The tube formulas in this paper may also be interpreted as fractal analogues of the 
Steiner formula and its generalizations (see Section 1.1 below, in particular (1.6)). Steiner- 
type formulas express the volume of the e-parallel sets of a given set A as a polynomial 
in e with coefficients that just depend on A. Under the additional assumption that the 



We denote by res (/(.v); .v = ui) the residue of a meromorphic function / at an isolated singularity u). Recall 
that this is the unique value a such that (/(i) - a)l(.s - a)) has an analytic antiderivative in a punctured disk 
[s \ Q < \s - Ci)\ < S\\ equivalently, the residue is the coefficient a_i in the Laurent expansion of /. 




(1.3) 
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complex dimensions (i.e., the poles of are simple, our tube formula can be written in 
the following way (see Corollary 5.9): 



where the coefficients q and et(e) are given (in (5.19) - (5.21)) in terms of the residues 
of ^(s) and the geometry of the generator G. If, in addition, the tube formula of the 
generator G is a polynomial, then the coefficients e^is) disappear (see Corollary 5.10) and 
the remaining coefficients are completely independent of s, just as ffie coefficients in the 
Steiner formula. (Compare the "fractal power series" in formula (5.23) to the polynomial 
in (1.6).) 

This paper is part of the program of the present authors to develop a fractal notion of 
curvature in terms of complex dimensions, and to relate it to other notions of curvature, 
especially as developed in [Win, LlWi]. 

We note that related questions are also being concurrently studied by other researchers 
[DKOU]. Recently, some tube formulas extending aspects of [LaPe2, LaPe3] have been 
obtained for tilings associated to graph-directed iterated function systems in [DDKU]. 

For our purposes, the precise embedding of T into is not important and the mapping 
*Py associated to is not emphasized. Due to the disjointness of the sets G^ in Defini- 
tion 1.2, the tube formulas require only those properties of fractal sprays which depend 
either on the geometry of the generator G or on the scaling ratios fj. 

In particular, for the generator G, we will require that the inner parallel volume of G 
admit a Steiner-like formula (Definition 2.1); that is, it can be represented as a 'polynomial' 
in e where the coefficients are allowed to depend on e. The Steiner-like condition should 
not be viewed as a restriction on the class of allowed generators G but as a choice of 
the representation of its inner parallel volume. In particular, Steiner-like representations 
are not unique. For the fractal string we assume that it is languid or strongly languid 
(Definition 3.3 or Definition 3.4), which is similar to the assumptions made in previous 
tube formula results. In the case of self-similar tilings, these languidity assumptions are 
always satisfied. We describe these conditions in detail in the following sections. 

Remark 1.3. Without loss of generality, and in contrast to [LaPe2, LaPe3, Pea, PW], we 
make a normalization assumption on the fractal string, for the remainder of the paper: 



This assumption imposes no restrictions on the class of fractal sprays, but will simplify 
the exposition greatly. It amounts to choosing the largest connected set in the spray as the 
generator (or one of them, if there is not a unique largest set). Also, instead of thinking 
of the numbers (j as distances (as in [LvF2], where the terms in a fractal string represent 
usually lengths of subintervals of R), we think of them as scales or scaling ratios. Thus, 
€\ is the scaling factor of the identity mapping I : R'' ^ R'', in accord with the original 
definition of fractal sprays given in [LaPo2] (see also [LvF2]) and the interpretation in 
terms of self-similar tilings discussed in Section 5 and in [LaPe3, PW, Pea]. 

1.1. Tube formulas and classical geometry. To motivate our theorems, we give a brief 
description of tube formulas in geometry. Such formulas have myriad applications in con- 
vex, integral and differential geometry and have roots in the results of Steiner [Ste] (when 
A is convex) and Weyl [Wey] (when A is a smooth submanifold). For connections to con- 
vex and integral geometry, see [KR, Sch2], and for connections to diflFerential geometry. 



d 




(1.4) 
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see [BG, Gra]. For a bounded set A c R'' and e > 0, we denote the s-neighborhood (or 
E-pamllel set) of A by 

As := {x € R"' \ A ■■. dist(x, A) < e). (1.5) 

Sometimes A^ is referred to as a "collar" in the hterature. Note that some authors include 
the set A in A^, but we have instead excluded A from Ag. In particular, A^ is not a neigh- 
borhood of A in the topological sense. 

The Steiner formula is a foundational result of convex geometry which states that the 
tube formula of any compact convex subset of R'' is a polynomial in s. 

Theorem 1.4 (Steiner formula). IfK c R"' is convex and compact, then the d-dimensional 
volume of Kb is given by 

d-\ 

UKe) = Yj ^'''(^d-kVkiK), (1.6) 

where the coefficients VkiK) depend only on the set K, and aj is the volume of the unit ball 
in W. 

Note that formula (1.6) can be extended to 

d 

Aa(K VJK,)^Yj s'^'a^^kYkiK), (1.7) 

k=0 

where Vd(K) :- Ad{G) is the volume of K. The coefficients Vo(K), . . . , Vd{K) are called in- 
trinsic volumes or Minkowski functionals of K. They form a system of important geometric 
invariants which is, in a sense, complete. Some of them have a simple direct interpretation. 
In particular, is the volume of K and Vd-\ is half the surface area of its boundary (pro- 
vided K has interior points); furthermore, Vi is, up to a normalization constant, the mean 
width of K, while Vq is its Euler characteristic. For nonempty convex sets K, Vo is always 
equal to 1. See [Sch2, Section 4.2] for further details. 

When the set K is sufficiently regular (i.e., when its boundary is a surface), the 
coefficients VtiK) can be given in terms of curvature tensors, and the Steiner formula 
coincides with the tube formula obtained by Weyl in [Wey] for smooth submanifolds of R'' 
without boundary. In [Fed], Federer unified both approaches and extended these results to 
sets of positive reach through the introduction of curvature measures Co{K, ■),..., Cd{K, ■) 
and a localization of the Steiner formula. A set K c M.'^ is said to have positive reach iff 
there is some 6 > such that any point x e W' with dist(jic, K) < 6 has a unique metric 
projection Pk(x) to K; i.e., there is a unique point pxix) in K minimizing dist{x, K). The 
supremum of all such numbers 5 is called the reach of K. 

The intrinsic volumes Vk{K) turn out to be the total masses of the curvature measures: 
VkiK) = Ck{K,W') fork = Q,...,d. Here, the volume measure Cd{K, •) := Ad{K n ■) is 
added for completeness. 

Federer's curvature measures and associated tube formulas have since been extended 
in various directions; see, for example, [Schl, Zahl, Zah2, Fu, RZl, RZ2] along with the 
book [Sch2] and the references therein. Recently (and most generally), so-called support 
measures have been introduced in [HLW] (based on results in [Sta]) for arbitrary closed 
subsets of R'', which were also shown to admit a Steiner-type formula. 

The total curvatures and the curvature measures above are defined as the coefficients of 
some tube formula. It is precisely this approach that we hope to emulate in a forthcoming 
paper, making use of the tube formulas obtained in the present paper We believe that 
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(for a suitable choice of the Steiner-Hke representation for the generators) the coefficients 
appearing in our tube formulas may also be understood in terms of curvature, in a suitable 
sense, and that a localization of the results in this paper may lead to a notion of complex 
curvature measures (or possibly, distributions) for fractal sets. We hope to explore such a 
possibility in a future work; see also Section 8.6. 

1.2. Outline. In Section 2, we discuss the geometric hypotheses placed upon the gener- 
ator(s) of the fractal spray. In Section 3, we define a zeta function associated to a fractal 
spray T, which will allow us to derive a pointwise tube formula for T in Section 4. In 
Section 5, we obtain the tube formula associated with a self-similar tiling (an important 
special case of a fractal spray). Several examples are discussed in Section 6. In Section 7, 
we give the detailed proof of the main theorem (Theorem 4.1), the pointwise tube for- 
mula for fractal sprays, as well as of Corollary 5.9, the fractal tube formula for self-similar 
tilings. Finally, in Section 8, we discuss the relation with previously obtained tube formulas 
and give some possible directions for future research. 

2. Steiner-like formulas for generators. 

In this paper, we consider the interior of a set instead of its exterior, as discussed in 
Section 1.1. However, our primary requirement of a generator is that it has a similar (in- 
ner) tube formula; see Definition 2.1 below and also Section 1.1 for motivation of the 
nomenclature. 

For a nonempty and bounded open set G c R'', let g - p{G) denote its inradius; that is, 
the radius of the largest open ball contained in G. It is clear that g is always positive and 
finite. In case G is the generator of a fractal spray 7~, we have 

p{G^) = p^'jG) = p({jG) = EjpiG) = Ejg (2.1) 

for the inradii of the components G^ of T. 

It will be useful to write the inner parallel volume V{G, s) of the set G c R'' as a 
"polynomial-like" expansion in e of degree at most d. More precisely, we have the follow- 
ing definition. 

Definition 2.1. An (inner) Steiner-like formula (or a Steiner-like representation of the 
tube formula) for a nonempty and bounded open set G c R'' with inradius g - p(G) is an 
expression for the volume of the inner e-parallel sets of G of the form 

d 

V{G, s) = Yj '^kiG, s)s''-'', for < s < g, (2.2) 

k=0 

where for each k - 0,1, . . . ,d, the coefficient function Kit(G, ■) is a real-valued function on 
(0, g] that is bounded on [sq, g] for every given sq e (0, g]. 

Remark 2.2 (The choice of the coefficient functions K^iG, e)). Note that a representation 
of the form (2.2) always exists. For example, one can always take a trivial representation 
with Kj{G, s) - V(G, s) and ko{G, s) - ■ ■ ■ = Kd-\{G, e) = on (0, g}. Another, slightly less 
trivial, representation is given by letting Kk{G,e) - ■j-^V{G,b)£!^^'' for ^ = Q,. . . ,d. For 
brevity, we may use the term "Steiner-like generator/set" to indicate that a fixed Steiner- 
like representation for the tube formula is intended, and write "tube formula" for "inner 
tube formula". 

We have in mind nontrivial representations of the volume function, in which the coef- 
ficients allow interpretations in terms of curvature. Clearly, not every representation can 
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have such an interpretation, and so some uniqueness condition will be needed to character- 
ize the correct one for this purpose. However, this is not our aim here (this issue shall be 
addressed in a forthcoming paper by the same authors). For the main results of this paper, 
the tube formulas for fractal sprays (and tilings), we make no assumptions on the genera- 
tors. In fact, our theorems provide many tube formulas for the same spray — one for each 
choice of a Steiner-like representation for the generators. Our formulas should be seen as 
a general tool to transfer a given representation of the volume function of a generator into 
a tube formula for the generated sprays. We do not yet know what the canonical choice of 
the representation for the generator is, but our approach seems flexible enough to contain 
it. It seems that a reasonable strategy would be to choose the coefficients as "constant as 
possible". It is likely that some integrals of the support measures of [HLW] could provide 
the coefficients of some canonical representation. 

Remark 2.3 (Monophase and pluriphase generators). As noted above, the coefficients in the 
expansion (2.2) are clearly not unique. However, if a set G has a Steiner-like representation 
with constant coefficients 

Kk(G, s) = KkiG) for all s e {0, g] and k = 0,1, . . . ,d, 

then such an expansion is unique, and the set is called monophase. More precisely, a 
bounded open set G c R'' is monophase iff its inner tube formula may be written in the 
form 

d 

V(G, s)^Yj '<k(G)e''-', 0<s<g. (2.3) 

For monophase sets, we always choose this canonical Steiner-like representation. In this 
case, one has KjiG) = 0, since otherwise limg^o ^^diG-s) + 0. The monophase case has 
been treated in [LaPe3], at least from the distributional perspective. A variety of natural 
and classical examples of self-similar tilings in R'' have monophase generators; see [LaPe3, 
Section 9]. Furthermore, all geometric (or ordinary) fractal strings (i.e., 1 -dimensional 
fractal sprays) also have monophase generators, since G is always an interval; see [LaPe3, 
Section 8.2]. In general, however, it is rather restrictive to assume that the generator is 
monophase because many sets (including generators of self-similar tilings) do not have a 
polynomial expansion with constant coefficients; see Section 6. 

For monophase sets, the inner tube formula is a polynomial for e e (0, g\ and this is 
the reason for the nomenclature. More generally, as in [LaPe3, LaPe2], we say a bounded 
open set G c R'' is pluriphase iff it has a Steiner-like tube formula with coefficient func- 
tions /Ci(G, ■) that are piecewise constant with respect to a finite partition of [0, g\. In short, 
the inner tube formula is piecewise polynomial, with finitely many pieces. (Such a repre- 
sentation is unique if it is assumed that one takes the partition to have as few components as 
possible.) We use the term general Steiner-like (or Steiner-like with variable coefficients^ 
to emphasize the distinction from the special cases of monophase and pluriphase sets. It 
was conjectured in [LaPe3, LaPe2] that all convex polyhedra are pluriphase. 

Remark 2.4. The above definition of "Steiner-like" is slightly more general than in [LaPe3, 
Definition 5.1], where it was introduced: in particular, the local integrability of each coeffi- 
cient function ki(G, ■ ) and the limit condition for lim£^o+ K-k{G, s) have both been removed. 
The content of Proposition 2.5 was taken as a hypothesis in [LaPe3], but is now seen to 
follow from the assumptions in Definition 2.1. 
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No assumption is made on the uniqueness of the coefficients Kk{G,e) in Definition 2.1 
(as discussed in Remark 2.2), but any choice of coefficients for G satisfying (2.2) gives rise 
to some coefficients kic(G\ e) for each set G^ = ^ jiG) by defining 

Kk(G^, e) := {''jKkiG, ifs), for < e < p(G-') = (jg, (2.4) 

as is seen in the following proposition. Here and henceforth, denotes the A;* power of (j. 

Proposition 2.5. Let T — {G-') be a fractal spray on a generator G with a given Steiner- 
like representation as in (2.2). Then the inner tube formula of each set G' has a Steiner-like 
representation in terms of the same coefficients: 

d 

V(Gj, e) = 2 {'jKkiG, {fs)e''-\ forQ<s< piG^) = {jg. (2.5) 

k=0 

Proof. The motion invariance and homogeneity of Lebesgue measure implies that for each 
j > 1, A^iGi^;) = Ad{^'j(G-l;/Ci)) - ijA^iG-E/fj), where ^'j is the similarity transformation 
of R'' described in Definition 1.2. Whence, by (2.2) and (2.4), 

d d 

ViGj, s) = £''jV(G, i-j's) ^{'jYj •'kiG, q'sW-j's)'-' = ^ ^,(G^ e)e''-^ 

k=Q k=0 

for < e < p{G^) = £jg, and the Steiner-like representation (2.5) follows. Note that the 
coefficients Kk{G^ , ■) of G^ clearly inherit the boundedness properties from the coefficients 
KkiG, ■) of G. □ 

In the sequel, we will always work with the coefficient functions of the sets G-' chosen 
according to (2.4). Proposition 2.5 ensures this choice is always possible. 

Up to this point, the coefficient functions Kk{G, ■) in a Steiner-like formula for G have 
been defined only for Q < e < g - p(G). For k - . . ,d, we define /<'jt(G) := Ki^{G,g) 
and then extend Kk(G, s) to e e (g, oo) as constant functions with this value: 

Kk(G, e) KkiG) for e>g. (2.6) 

Note that (2.2) need not hold for e > ^ and so we have the freedom of the choice (2.6). We 
emphasize that this choice is vitally important for the tube formulas in Theorem 4. 1 and its 
corollaries below to be correct. 

Note that, for e - g, equation (2.2) implies that the i/-dimensional Lebesgue measure of 
G satisfies 

d d 

AAG) = y(G, g)^Yj S'^S"'' = Z ^*(^)^'"'- (2.7) 

*=0 k=0 
3. ZeTA functions and COMPLEX DIMENSIONS 

We will require certain mild hypotheses on the fractal string X which gives the scaling 
of the spray T. These conditions are phrased as growth conditions on a zeta function 
associated with within a suitable window, as defined just below. 

Definition 3.1. For a fractal string £. - {(jfj^^^, the scaling zeta function is given by 

oo 

= (3.1) 

for i e C with Re(i) > D, where D is the abscissa of convergence of this series. (Compare 
with [LvF2, Def. 1.8], where ^ is called the geometric zeta function of £. ) Recall that 
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D :- inf{a e R . TjJ^i < ^nd that ^ is holomorphic (i.e., analytic) for Re s > D. 
Henceforth, if W c C contains {Re s > D] and ^ has a meromorphic continuation (nec- 
essarily unique) to a connected open neighborhood of W, we abuse notation and continue 
to denote by its meromorphic extension. Under these assumptions, each pole oj e W of 
^ is called a visible complex dimension of X, and the set of visible complex dimensions is 
written as 

D_c(W) = {w € w is a pole of ^^). (3.2) 

Moreover, in the special case when ^ has a meromorphic continuation to all of C, we 
may choose W - C and then simply write := T)_£(C) and refer to as the complex 
dimensions of £, . 

In practice, W will be a window (the part of C to the right of a screen 5 ) as in Definition 
3.2, just below. The following three definitions are excerpted from [LvF2, Section 5.3]. 

Definition 3.2. Let 5 : R ^ D] be a bounded Lipschitz continuous function. Then 
the screen is 5 = {^(f) + if •. t e R), the graph of a function with the axes interchanged. 
Here and henceforth, we denote the imaginary unit by i := V-T. We let 

mfS := inf,<=R5(0 = inf{Re s ■. s e S), and (3.3) 
sup 5 :- sup,gjj 5(f) - sup{Re s -.s&S). (3.4) 

The screen is thus a vertical contour in C. The region to the right of the screen is the set 
W, called the window. 

W:={zeC-. Rez>5(Imz)). (3.5) 

For a given string we always choose S to avoid 2)x ™d such that ^ can be meromor- 
phically continued to an open neighborhood of W. We also assume that sup S < D, that is, 
S(t) < D for every f e R. 

Definition 3.3. The fractal string X. is said to be languid if its associated zeta function 
^ satisfies certain horizontal and vertical growth conditions. Specifically, let {r„)„gz; be a 
sequence in R such that r_„ <0 <T„ for n > 1, and 

T 

lim T„ = oo, lim r_„ = -oo, and lim ■ — ^ = 1. (3.6) 

For X to be languid, there must exist constants y e R and c > 0, and a sequence {T„} as 
described in (3.6), such that: 

For all n e Z and all a > S(T„), 

\(^(a + iT„)\<c-{\T„\ + iy, LI 

and for all f e R, |f| > 1, 

|^^(5(f) + if)|<c-|fr. L2 
In this case, X. is said to be languid of order j . 

Definition 3.4. The fractal string X is said to be strongly languid of order y and with 
constant A iff it satisfies LI and the following condition L2', which is clearly stronger than 
L2: 

There exists a sequence of screens S,„ for m > 1, f 6 R, with sup5,„ — > -oo as m — > oo, 
and with a uniform Lipschitz bound, for which there exist constants y € R and A,c > Q 
such that 
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\^^(S„.(t) + it)\<c-A^''"^'H\t\ + iy, L2' 

for all f e R and m > 1. 

By saying is languid", we mean just that X. is languid. In the rest of this paper, 
7~ is assumed to be a fractal spray with a generator G c R'', scaled by the fractal string 
-C - {^j)J=i with = 1 as in Remark 1.3. The tubular zeta function first appeared in 
[LaPe3], but we need to modify the definition in order to extend it to the case when the 
generators are not monophase; thus, the following definition is new. 

Definition 3.5. The tubular zeta function of the fractal spray T is defined by 



j=i \k=a 



(3.7) 



for every s e (0, oo) and for each s e C such that the sum converges absolutely. As in 
Definition 3.1, we will henceforth abuse notation and use ^^-(e, s) to mean a meromorphic 
extension of the function defined by the formula (3.7), as convenient. 

Note that by (2.7), (3.1), and Proposition 2.5, for s > g, one has 



L^Ti 7-^ ^x(^)z. 



A=0 



d-l 



g'-'Kt(G){d-k) 



(3.8) 



Since ^jis) is a Dirichlet series, it has an abscissa of convergence: there is a unique number 
D e [-00, oo] such that ^(i) converges absolutely for s with Re i > D and diverges 
for s with Re s < D. The abscissa of convergence is thus analogous to the radius of 
convergence of a power series. Note that in all reasonable situations we have Q < D < d. 
Indeed, D > follows immediately from the non-finiteness of the fractal string (assumed in 
Definition 1.1), and D < d follows if one requires that the generated fractal spray T have 
finite total volume. (Note that for fractal sprays with infinite total volume, the question 
for a (global) tube formula does not make sense.) For the scaling zeta functions of the 
self-similar tilings discussed in Section 5, one has Q < D < d, and D coincides with the 
Minkowski dimension, Hausdorff dimension, and similarity dimension of the associated 
self-similar set; for a precise statement, please see [LaPe3, Section 4.3] and Section 5 
below (especially Remark 5.4). 

Note that the tubular zeta function may be viewed as a generating function for the 
geometry of the fractal spray T. 

Proposition 3.6 clarifies the relation between the scaling zeta function ^ and the tubular 
zeta function (j- of a fractal spray T. It also motivates and justifies the definition of com- 
plex dimensions of fractal sprays. The intended application of Proposition 3.6 is with Q. as 
a suitable open neighborhood of a window W for the scaling zeta function ^ of T, as in 
Definition 3.2. (Proposition 3.6 is extended significantly in Theorem 7.2 and Lemma 7.9.) 

Proposition 3.6. IfD is the abscissa of convergence of ^j-, then for all e > 0, the series in 
(3.7) defining ^■j-(_b, s) converges absolutely for any fixed s E C\[0,\, . . . ,d] with Re s > D. 
More generally, suppose ^ is meromorphic in a connected open set Q containing {Re s > 
D]. Then for all e > 0, the function ^j-ie, ■ ) is meromorphic in Q. and each pole to € Q, of 
t^<j-(e, ■ ) is a pole of^^ or belongs to the set {0,1, ... , d}. 
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Proof. Fix e > 0. Upon expanding (3.7) of Definition 3.5 and interchanging the sums, the 
tubular zeta function becomes 

(^(e, ^'M*^^ ^/^^ - ^j^(j:(s)Ad(G). (3.9) 

k=0 j=l 

It is clear that the second term on the right-hand side of (3.9) is convergent for s as in the 
hypotheses, so it remains to check that the first term is similarly convergent for each k. 
Since e is fixed and Hj \ 0, define J - 7(e) to be the index of the last scale greater than 

e: 

7(e) := max{; > 1 •■ (J^s <g]wO. (3.10) 

At the end of (3.10), "VO" indicates that /(e) = for e > tig? Now KkiCq^e) = KkiG) 
for all j > J, by (2.6), and so 

oo 7 oo 

;=i j=i j=J+i 

Observe that the first sum on the right-hand side of (3.11) is entire, as a finite sum of the 
entire functions Cj.x\ and that the second sum on the right-hand side of (3.11) converges 
absolutely exactly where ^ does; that is, for Re s > D. 

Justification of the claims of meromorphicity are obtained by parallel reasoning; the 
decomposition (3.11) shows that, except possibly for {0, 1, ... , d}, the tubular zeta function 
is meromorphic precisely where the scaling zeta function is. □ 

Definition 3.7 (Complex dimensions). The set 

= ©i:U{0, l,...,t/) 

of (potential) poles of ^j- is called the set of complex dimensions of T. Let W c C be a 
window for ^ as in Definition 3.2, so that ^ is meromorphic in an open neighborhood of 
W. (Proposition 3.6 thus implies that for each fixed e > 0, the function (ri^, ■ ) is also 
meromorphic in an open neighborhood of W.) Then Dt{W) := Dt n W is called the set of 
visible complex dimensions of X in W, in parallel with (3.2). We refer to as the scaling 
complex dimensions and {0, 1, . . . , as the integer complex dimensions of T. 

4. POINTWISE TUBE FORMULAS FOR FRACTAL SPRAYS 

Now we are ready to state one of our main results, a pointwise tube formula for a 
fractal spray T, which, for e > 0, describes the inner parallel volume V{T, e) as a sum 
of the residues of its tubular zeta function ^^-(e, s). For fractal sprays with more than one 
generator, one can consider each generator independently, and the tube formula of the 
whole spray is then given by the sum of the expressions derived for the sprays on each 
single generator Thus, there is no loss of generality in considering only the case of a 
single generator in Theorem 4. 1 . 

Theorem 4.1 (Pointwise tube formula for fractal sprays). Let T be a fractal spray given by 
the fractal string X = '^^d the generator G C R''. Fix a Steiner-like representation 

for G, as in (2.2), and assume that the abscissa of convergence D of the scaling zeta 
function ^ ofT is strictly smaller than d. 



By convention, sup = -oo. Thus the maximum VO is included so that /(e) = when e > tig. 
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Tube formula with error term. If ^ is languid of order y < \ for some screen S for 
which 5 (0) <Q (so that W contains the integers [0,1, ... ,d]}, then for all e > 0, 

y(r,£)= Yj res(^^fei);i = + + (4.1) 
where the error term "R (given explicitly in Remark 4.3 below) is estimated by "Ris) — 

Tube formula without error term. If ^ is strongly languid of order -y < 2 and with 
constant A > 0, then the choice W — C for the window is possible in (4.1), implying that 
the error term 'R(s) vanishes identically for all < e < mm{g,A^^g}. 

Theorem 4. 1 and its corollaries are consistent with earlier results in [LvF2] and [LaPe3] 
and generalize them in several respects; see Section 8.1. We will give the rather lengthy 
proof of Theorem 4.1 in Section 7.1. For a description of one of the main new ideas and 
techniques, we refer the reader to Remark 7.1. 

Note that the tube formula without error term is an exact pointwise formula. For this 
result, one must assume that the sequence of screens of Definition 3.4 satisfies 

Sm(0) < for each m. However, there is no loss of generality because sup5,„ — » -oo. 

The following result is really a corollary of the proof of the first part of Theorem 4.1. 
For this reason, its short proof is provided in Section7.1.9 atthe very end of Section 7.1. 

Corollary 4.2 (The monophase case). If in addition to the hypotheses of the first part 
of Theorem 4.1, we assume that G is monophase, then the tube formula with error term 
remains valid ( with the same error estimate ), without the restriction that S (0) < 0, provided 
this hypothesis is replaced by the much weaker condition that the screen S avoids the 
integers 0,1, ... ,d. Hence, in particular, it still holds for a screen S that is arbitrarily 
close to the vertical line Re s — D. Moreover, the error term 'R is given by (4.3) (or (4.2)j. 

Remark 4.3. The error term H in formula (4. 1) in Theorem 4. 1 is explicitly given by 



2ni Js d-s f-^ 



\k=o ^ ■ 



-(d - k)Ku(G) 



ds. (4.2) 



The integrand in formula (4.2) will be called the tail zeta function of T and denoted by 
^7"iaii(^' ™ Section 7; see, in particular. Section 7.1.1 and equation (7.4). The function 
^rtaii(^' *) is oii^ of the head-tail spUtting of the tubular zeta function i^j-ie, s) employed 
in the proof of Theorem 4. 1. In the situation of Corollary 4.2, one has = i^j- and thus 
the error term is equivalently given by 

Kis)^ ^^^ris,s)ds. (4.3) 

See Section 8.1 below for a discussion of the consistency of the error term with earlier 
results. 

Remark 4.4. For investigating delicate questions concerning the Minkowski measurabil- 
ity of fractal sprays and self-similar tilings (see, for example, [LaPe3, Corollary 8.5]), it 
is important to be able to drop the assumption that 5(0) < 0, as in Corollary 4.2. How- 
ever, this generalization poses technical challenges for the case of more general (i.e. non- 
monophase) generators. In the monophase case, in contrast, our tube formulas enable us 
to derive results on the Minkowski measurability of fractal sprays. For example, for a 
self-similar tiling T (as discussed in Section 5 below), let us denote 



yt=0 
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Assume that ToiG) + and (in the lattice case) that ro+i,„p(G) + for some m eZ \ {0}? 
If d - I < D < d, and the self-similar tiling T has a single monophase generator, then 
one can apply the methods of proof (and the conclusions) of Theorems 8.23, 8.30, and 
8.36 (along with Theorems 2.17, 3.6, 3.25 and 5.17) of [LvF2] to see that T is Minkowski 
measurable if and only if ^ is nonlattice. In addition, 7~ has Minkowski dimension D. In 
the nonlattice case, T has (positive and finite) Minkowski content 

resfe(.);D)E^^, (4.5) 

with residues res (^(i); Dj given by (5. 15). In the lattice case, T has average Minkowski 

content given by (4.5), with residues res (^(i); £>) as in (5. 14). 

With definitions suitably adapted from [LvF2, Chapter 8], an entirely analogous state- 
ment can be made about a self-similar set F. More precisely, if the tiling is also assumed 
to satisfy the compatibility condition (5.7), then analogous results extend to the associated 
self-similar fractal F. This strengthens and specifies the results of [LaPe3, Corollary 8.5]; 
see also [LaPe3, Remark 10.6]. Further discussion of this issue is lengthy and beyond the 
scope of the present paper. For the interested reader, details on the monophase case can 
be found in [LPW2]; see also Section 8.4. For more general generators, these results are 
under further development in [LPWl]. 

Remark 4.5. In light of Definition 3.3 and Definition 3.4, note that if ^ is strongly languid 
of order y, then it is also strongly languid (and hence languid) of any higher order, but 
not necessarily of any lower order Consequently, the assumptions of the second part of 
Theorem 4. 1 (the strongly languid case) do not imply those of the first part (the languid 
case). Compare to [LvF2, Remark 8.8]. 

Remark A.6. Define T :- IJ^i and let r_£ be as defined in (1.1). Note that homogeneity 
of Lebesgue measure gives 

oo oo 

and hence the tube formula (4.1) expresses the fact that the measure of the complement of 
r_e in T is given by 

^^(r\r_,) = - Yj res(^^(e,i);i = w) ( - 9?(e)). (4.6) 

5. POINTWISE TUBE FORMULAS FOR SELF-SIMILAR TILINGS 

In [Pea, PW, LaPe3], the focus is on self-similar tilings. Such an object is a fractal 
spray associated to an iterated function system {Q>\, . . . ,(^n], N > 2, where each O,, is 
a contractive similarity mapping of with scaling ratio r„ e (0, 1). For A c R"', we 
write 0(A) := U^^=i '^n(A)- The self-similar set F (generated by the self-similar system 
{Oi, . . . ,(!>ff}) is the unique (compact and nonempty) solution of the fixed-point equation 
F = 0(F) ; cf. [Hut]. The fractal F is also called the attractor of the self-similar system 
{Oi, . . . ,^m]- To proceed with the construction of a self-similar tiling, the system must 
satisfy the open set condition and a nontriviality condition: 

A self-similar system (Oi, . . . ,0^r) (or its attractor F) satisfies the open set condition 
(OSC) if and only if there is a nonempty open set O QW' such that 



'See Remark 5.6 for a discussion of terminology and notation in the lattice and nonlattice cases. 
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(t>„(0)cO, n^l,2,...,N (5.1) 
<D„(0) n (D,„(C>) = for « m. (5.2) 

In this case, O is called a. feasible open set for {<l>i, . . . , (Da;) (or F), cf. [Hut, Fal, BHR]. 

A self-similar set F satisfying OSC is said to be nontrivial, if there exists a feasible 
open set O such that 

Ot^iO), (5.3) 
where O denotes the closure of O; otherwise, F is called trivial. This condition is needed 
to ensure that the set O \ 0((9) in Definition 5.2 is nonempty. It turns out that nontriviality 
is independent of the particular choice of the set O. It is shown in [PW] that F is trivial 
if and only if it has interior points, which amounts to the following characterization of 
nontriviality: 

Proposition 5.1 ([PW, Corollary 5.4]). Let F Q be a self-similar set satisfying OSC. 
Then F is nontrivial if and only if F has Minkowski dimension (or equivalently, Hausdorff 
dimension ) strictly less than d. 

All the self-similar sets F considered in this paper will be assumed to be nontrivial, and 
the discussion of a self-similar tiling T implicitly assumes that the corresponding attractor 
F is nontrivial (and satisfies OSC). 

Denote the set of all finite words formed by the alphabet {1, ... ,A^) by 

oo 

W |J{1,...,A')*. (5.4) 

For any word w = wiW2 ...Wn e "W, let r„ :- r,,,, ■ . . . ■ r,,,, and O^. :- O^.j o ■ • ■ o In 
particular, if w e W is the empty word, then r„, = 1 and - Id. 

Definition 5.2. (Self-similar tiling) Let (9 be a feasible open set for {Oi , . . . , Oat). Denote 
the connected components of the open set O \ <1>((5) by Gq, q & Q. Then the self-similar 
tiling T associated with the self-similar system {Oi, . . . , 0^) and O is the set 

T(0) := {0,,(G,) -.we^^qe Q}. (5.5) 

We order the words w^^\ w^^\ ... of 'W in such a way that the sequence {{j}J^i given 
by £j :- r„u), j - 1,2, . . ., is nonincreasing. It is clear that a self-similar tiling is thus a 
collection of fractal sprays, each with fractal string £. - {£j]°J^y and a generator Gq, q ^ Q. 
In this context, the mapping 'V j appearing in Definition 1.2 corresponds to O^oj. 

The terminology "self-similar tiling" comes from the fact (proved in [PW, Theorem 5.7]) 
that T{0) is an open tiling of O in the following sense: The tiles <i>H.(Gg) in T{0) are pair- 
wise disjoint open sets and the closure of their union is the closure of O, that is. 

Remark 5.3 (Tube formulas for self-similar sets). In [PW, Theorem 6.2], precise conditions 
are given for when the tube formula of a self-similar tiling can be used to obtain the tube 
formula for the corresponding self-similar set, the attractor F; recall from (1.5) that for a 
bounded set A c R'', we define :- {x e M.'' \ A ■. dist(jc, A) < e}. 

Let F c R'' be a self-similar set satisfying OSC with some feasible open set O and 
diniM F < d (i.e., F is nontrivial). Let T(0) be the associated tiling of O, and let K :- 
O and r := U 

we'W.qeQ ^w(Gq). Then [PW, Theorem 6.2] states that one has a disjoint 

decomposition 
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= U K^, for all e > 0, (5.6) 
if and only if the following compatibility condition is satisfied; 

hAKQF. (5.7) 



In this case, the tube formula for the self-similar set F can be obtained simply by adding 
to V{T,e) the (outer) tube formula AdiKe) as in (1.6) (although note that in the present 
context, K need not be convex). For example, the Sierpinski gasket and the Sierpinski 
carpet tilings (see Figures 6.2 and 6.3) satisfy the compatibility condition (5.7), whereas 
the Koch curve and the pentagasket tilings do not (see Figures 6.1 and 6.4). Condition 
(5.7) will not be assumed in the remainder of the paper. 

From now on, let T - T(0) be a self-similar tiling associated with the self-similar 
system {'^n}^^i and the generator G. We refer to the fractal F as the self-similar set as- 
sociated to T. For the same reasons as described in the first paragraph of Section 4, we 
lose no generality by stating all results for self-similar tilings with one generator, which we 
will denote by G in the sequel. (For natural examples of a self-similar tiling with multiple 
generators, see the pentagasket depicted in Figure 6.4 of the examples section, and also 
Example 6.2, which is depicted in Figure 6.8.) Without loss of generality, we may also 
assume that the scaling ratios {'"nl^^i of {<I>„)^^j are indexed in descending order, so that 

< r/v < ■ ■ ■ < r2 < ri < 1. (5.8) 

It follows from [LvF2, Theorem 2.9] (see also [LaPe3, Theorem 4.7]) that ^ has a 
meromorphic extension to all of C given by 

= — ^eC, (5.9) 

and hence that the set of scaling complex dimensions of T consists precisely of the 
roots 5 e C of the equation 

N 
n=I 

It is known from [LvF2, Theorem 3.6] that the set lies in a bounded vertical strip: there 
exists a real number D; e (-oo, D) such that 

Di < Re .? < D, for all s e D^. (5.11) 
For the remainder of this paper, we let 

Dr ^Dr(C)^D_c^ {0,1, ■■■,<!}. (5.12) 

Remark 5.4 (Various incarnations of D). Recall that D denotes the abscissa of convergence 
of It follows from [LvF2, Theorem 3.6] that D - D£ is a simple pole of ^ and that 
D is the only pole of ^ (i.e., the only scaling complex dimension of T) which lies on 
the positive real axis. Furthermore, it coincides with the unique real solution of (5.10), 
often called the similarity dimension of F and denoted by 6. Since F satisfies OSC, D 
also coincides with the Minkowski and Hausdorff dimension of F, denoted hy Df and 
Hf, respectively. (For this last statement, see [Hut], as described in [Fal, Theorem 9.3].) 
Moreover, it is clear that D > since N > 2, and that D < d; in fact. Proposition 5.1 
implies D < d. In summary, we have 



16 



MICHEL L. LAPIDUS, ERIN P. I. PEARSE, AND STEFFEN WINTER 



0<D<d and D ^ S ^ Dp = Hp. (5.13) 

The following result is an immediate consequence of [LvF2, Theorem 3.6], which pro- 
vides the structure of the complex dimensions of self-similar fractal strings (even for the 
case when D may be larger than 1). 

Proposition 5.5 (Lattice/nonlattice dichotomy, see [LaPe3, Section 4.3]). 

Lattice case. When the logarithms of the scaling ratios r„ are each an integer multiple 
of some common positive real number, the scaling complex dimensions lie periodically on 
finitely many vertical lines, including the line Re i = D. In this case, there are infinitely 
many complex dimensions with real part D. 

Nonlattice case. Otherwise, the scaling complex dimensions are quasiperiodically dis- 
tributed (as described in [LvF2, Chapter 3]j and s — D is the only complex dimension with 
real part D. However, there exists an infinite sequence of simple scaling complex dimen- 
sions approaching the line Re i = D from the left. In this case (cf. [LvF2, Section 3.7. l]j, 
the set {Re s-.se D} appears to be dense infinitely many compact subintervals of [Di, D], 
where Di is as in (5.11). 

Remark 5.6. It follows from [LvF2, Theorem 3.6] that in the lattice case (i.e., when r„ - 
r''", n = 1, . . . , A^, for some < r < 1 and positive integers {A;„)^^j), the scaling complex 
dimensions have the same multiplicity and a Laurent expansion with the same principal 
part on each vertical line along which they appear. In particular, since D is simple (see 
Remark 5.4), all the scaling complex dimensions {D + impl^ez (where p - In/ log r^^) 
along the vertical line Re s = D are simple and have residue equal to 

res(4>(i);£)) = ^ . (5.14) 

In the nonlattice case, D is simple with residue 

rcs{(^(s)-D)^-j^ -. (5.15) 

Note that (5.15) is also valid in the lattice case. Proposition 5.5 and the contents of this re- 
mark are used when applying Theorem 5.7 and Corollary 5.9 to the examples in Section 6. 

5.1. Exact pointwise tube formulas. The following result is a consequence of the strongly 
languid case of Theorem 4. 1 when applied to self-similar tilings. 

Theorem 5.7 (Exact pointwise tube formula for self-similar tilings). Assume that T is a 
self-similar tiling with generator G C Mf^, and that a Steiner-like representation has been 
chosen for G as in (2.2). Then for all s € (0, g), 

V(T, e) = ^ res (^^(e, s); s ^ co) + AAGX^rid). (5.16) 

Proof. The open set condition and nontriviality condition ensure that D < d. It remains 
to show that ^ is strongly languid of some order y < 2 with constant A - r^. Indeed, in 
view of (5.8) and (5.9), ^ is strongly languid of order y = < 2 with constants A - rt^ 
and C = 1 > 0, as in Definition 3.4: 

= ,_L s ^ ('•^T'""" ' ' (5.17) 
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1 Clearly, it follows from (5.11) that the sequence of screens {S,„}'^,^i in Definition 3.4 may 

2 be chosen to be a sequence of vertical lines lying strictly to the left of min{Df,0) and 

3 tending to -oo. In particular, this ensures sup^,,, < for all m - 1,2, Applying the 

4 second part of Theorem 4. 1 with A - r^, we deduce that the tube formula for T has no 

5 error term and is given by (5.16) for all < e < min{g, r^^g} = g- (Note that since < 1, 

6 we have r^^g > g.) a 

7 Remark 5.8. Theorem 5.7 generalizes to higher dimensions the pointwise tube formula for 

8 self-similar strings (i.e. 1-dimensional self-similar tilings) obtained in [LvF2, Section 8.4]. 

9 The formula (5.16) holds pointwise, as opposed to the corresponding result in [LaPe3, 

10 Theorem 8.3], which was shown to hold only distributionally, also generalizes the tube 

11 formula for self-similar tilings (obtained in [LaPe3, Theorem 8.3]) to generators which 

12 may not be monophase (or even pluriphase; see Remark 2.3). 

13 Concerning the proof of Theorem 5.7, see also the discussion in [LvF2, Section 6.4] 

14 regarding the self-similar string X. = {ij]J=i (this is a generalized self-similar fractal string, 

15 in the sense of [LvF2, Chapter 3]). The following more explicit form of Theorem 5.7 is 

16 used to compute examples in Section 6. 

17 Corollary 5.9 (Fractal tube formula). Assume, in addition to the hypotheses of Theo- 

18 rem 5.7, that the poles of the tubular zeta function C/j- are simple {which implies that I>£ 

19 and {0, 1, . . . , c/) are disjoint). Then for all < e < g, we have the following exact tube 

20 formula: 

d 

V{T, e) = 2 c^-" + 2 (q + euis)) e''-\ (5.18) 

resfc(.v);.v^c^) ^ ^-Ar(j_fe) 

Co; := -, — Kk{G), forueDz, (5.19) 

d — CO j—' cj - k 

A— 

Q := Kt(G)^j.(k), for ke {0,1,..., d}, (5.20) 

ek(s):^Y,^){Kt(G,{]'e)-Kt(G)), for k e {0,1, ... ,d}, (5.21) 

22 and 7(e) := maxfy > 1 ■. ^y'g < g] V as in (3.10). Alternatively, one has 

d 

ViT, e) = ^ c„e''-'- + ^ Ck(s)e''-\ (5.22) 

ojeDj: k=Q 

23 where c^ is as in (5.19) and Ck(e) :— q + eit(e) with q and ekie) as in (5.20)-(5.21), for 

24 ke{0,l,...,d}. 

25 The proof of Corollary 5.9 is postponed to Section 7.2, as it is technical and depends 

26 on the terminology and technique developed in the first part of Section 7 (the proof of the 

27 tube formula for fractal sprays. Theorem 4.1). Corollary 5.9 also allows us to recover the 

28 pointwise version of [LaPe3, Corollary 8.7], where the generator G was assumed to be 

29 monophase. 

30 Corollary 5.10 (Fractal tube formula, monophase case). In addition to the hypotheses of 

31 Corollary 5. 9, assume that G is monophase. Then, for all < e < g, we have the pointwise 

32 tube formula for self-similar tilings: 



21 where 
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d 

where (for co e and q (for k — 0,1, . . . ,d) are as in (5.19) and (5.20), respectively. 

Proof. When G is monophase, each function Kj^iG, ■ ) is constant and equal to Kk{G), and 
hence ek{e) = for each e > and k = 0,1, ... ,d, where eit(s) is as in (5.21). Conse- 
quently, Corollary 5.10 follows immediately from Corollary 5.9. □ 

Remark 5.11. For an arbitrary fractal spray T satisfying the hypotheses of the strongly 
languid case of Theorem 4.1, it follows from Theorem 4.1 (instead of Theorem 5.7), that 
(5.16) holds pointwise for all < e < mm{g,A^^g}. If, in addition, all the complex 
dimensions of T are simple, then one can deduce from Lemma 7.9 (as in the proof of 
Corollary 5.9) that (5.18) holds; see also Remark 7.10 in this regard. Moreover, if G is 
assumed to be monophase, then (5.18) takes the simpler form (5.23). A parallel remark 
holds (under the assumptions of the languid case of Theorem 4.1) for the tube formulas 
with error term considered in Section 5.2. 

5.2. Pointwise tube formulas with error term. 

Theorem 5.12 (Pointwise tube formula with error term for self-similar tilings). Assume 
that T is a self -similar tiling with generator G, and that a Steiner-like representation for 
G has been chosen. Let S be a screen such that S(0) < (so that all integer dimensions 
are visible) and let W be the associated window. Then for all e > 0, 

y(r,e)= Yj res (^^(fi,i);i = w) + ^<i(G)<£(flf) + ??(£), (5.24) 

where the error term 9?(e) is given explicitly as in (4.2) and satisfies "Ris) — Ois'^^^"^^), as 
0+. 

Moreover, if G is monophase, then this same conclusion holds without the assumption 
that 5(0) < 0, as long as S avoids the set {0, 1, . . . ,d).^ In addition, 'R(s) is equivalently 
given by (4.3) in this case. 

Proof. This follows immediately from the first part of Theorem 4.1, since the proof of 
Theorem 5.7 implies ^ is languid of order y = < 1 along any screen S . When G 
is monophase, the latter claim follows from Corollary 4.2. Finally, it follows from the 
second part of Remark 4.3 that in the monophase case, 'Ris) is equivalently given by (4.2) 
or (4.3). □ 

The following result is the exact counterpart of Corollary 5.9 (or of Corollary 5.10 when 
G is monophase). 

Corollary 5.13 (Fractal tube formula with error term). Assume, in addition to the hypothe- 
ses of Theorem 5.12, that the visible poles of the tubular zeta function are simple (which 
implies that I>£{W) and [0,1, ... ,d] are disjoint). Then for all e > 0, 

V{T, e) = 2 c^-'^ + Y^(ck + ek{e))8''-'' + Hie), (5.25) 
weDziW) te{0,l,...,£/)nw 

where the error term "Ris) is as in (4.2) and c^, q, Ck are as in (5. 19)-(5.21). 



In particular, this allows for a screen S which lies arbitrarily close to the vertical line Re s = D. 
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Figure 6. 1. The Koch curve tiling. 



Moreover, if G is assumed to be monophase, then (5.25) holds for any screen which 
avoids the set [0,1, ... , d), and the formula takes the simpler form 

y(r,e)= 2c^/-"+ eke''-' + Hie), (5.26) 

we©^(W) te|0,l,...,£/)niv 

with the error term 'R(b) as in (4.3). 

Proof This follows from Theorem 5. 12 by the same methods as in Corollary 5.9 (or Corol- 
lary 5.10, when G is monophase). □ 

Remark 5.14. The significance of the assumption S (0) < 0, and more importantly, the need 
for being able to omit it, is discussed in Remark 4.4 and Section 8.4. See also [LPWl]. 



6. Examples 

Firstly, it should be noted that Theorem 5.7 implies that all tube formula results for the 
examples of self-similar tilings of [LaPe3, Pea, PW] are now known to hold pointwise. 
This includes the Koch tiling (Figure 6. 1 and [Pea, Fig. 2 & 3]), the Sierpinski gasket tiling 
(Figure 6.2 and [Pea, Fig. 6]), the Sierpinski carpet tiling (Figure 6.3 and [Pea, Fig. 7]), 
the pentagasket tiling (Figure 6.4 and [LaPe3, Fig. 5]), the Menger tiling (Figure 6.5 and 
[Pea, Fig. 8]), and the three U-shaped examples from [PW, Fig. 3] (see Figure 6.8 for one 
of them). The tube formulas of the first three of these examples can be found in [LaPe3, 
Section 9]. 

In Figures 6.1-6.6 as well as in Figure 6.8, the following sets are shown from left to 
right. The set O is the initial open set of the tiling construction. (In all examples except 
the U-shaped one in Figure 6.8, O is the interior of the convex hull of the underlying self- 
similar set.) The second set shows the generator(s) of the tiling (or, more precisely, the set 
O \ 0((9)), while the subsequent ones give the first iterates of the generator(s) under the set 
mapping O. The right-most set always shows the union of all tiles of the generated tiling 
T(0). 

Of the self-similar tilings mentioned just above, only the Cantor carpet tiling and the 
U-shaped tiling will be studied in more detail below. Apart from illustrating how the tube 
formulas are applied in general, these two examples exhibit some important new features of 
the results obtained in this paper. Indeed, the Cantor carpet tiling (discussed in Section 6.1) 
has a generator which is not monophase (and not even pluriphase), a situation not covered 
by previous results. Furthermore, the U-shaped example (discussed in Section 6.2) has a 
generator which is itself fractal, in the sense that it has arbitrary small features and exhibits 
some kind of self-similarity. Finally, the binary trees discussed in Section 6.3 and the 
Apollonian packings discussed in Section 6.4 are natural examples of fractal sprays which 
are not self-similar tilings. 

6.1. The Cantor carpet tiling. We consider the self-similar tiling associated to the Carte- 
sian product C X C c r2 of the ternary Cantor set C with itself; see Figure 6.6. By abuse of 
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Figure 6.5. The Menger sponge tiling has a Steiner-like generator which 
is neither convex nor pluriphase; see the computations for the Cantor 
carpet in Section 6.1, for which the Menger sponge is a 3-dimensional 
analogue. 



notation, we denote the associated self-similar tiling by . The fractal C x C is constructed 
via the self-similar system defined by the four maps 

^j{x)^\x+lpj, ;=1,...,4, 

with common scaling ratio = 5, and points pj being the vertices of a square, as seen in 
Figure 6.6. Consequently, the corresponding string £.^2 - [ijYJ^y has scales 

3-[log4 3;]^ ;=1,2,..., (6.1) 

where {x\ is the integer part of x. 

The Cantor carpet tiling is discussed here because it has a generator G which is not 
monophase (and not even pluriphase), as seen in Figure 6.7 and formula (6.2). The inradius 
of the generator is g = p(G) - ij (3 V2), where C is the side length of the initial square (we 
set ^ = 1 in the sequel), and the relevant partition of the e-interval (0, g\ is 
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Figure 6.6. The Cantor carpet tiling C 



r 



= e,, < e < £i 



El < e < £2 



Figure 6.7. The generator of the tiling is not pluriphase. 



The tube formula for the generator of this tiling is given by the following Steiner-Uke 
(but clearly not pluriphase) representation: 



V(G,e) = 



(;r-8)e2 + 12V2ge, 



< £ < 



TIE - 4 arccos 

2 



V2' 



V2 

e > 



<£<. 



(6.2) 



Here, for si < e < £2, the constant term 8^^ = | in (6.2) gives the area of the four 
"protrusions" of G which are completely contained in G_£. By (6.2), we can take the 
coefficient functions K]S(j, s) to be 



K^){G, s) = 

Ki(G,e) = 
KiiG, e) = 




< e < 



V2' 



V2 



< £ < 



fl2V2g, 



2g 



(0, 



V2e2 - g\ 



0<e< 4, 
V2' 



V2 

< e < 4, 

V2' 

^ < « < ^, 



<£<, 



(6.3) 



(6.4) 



Since g - V2/6 and KkiG, g) - /ct(G) for A; = 0, 1 , 2, it follows that 

^o(G) = 0, Ki(G)^2g^^, and ^2(G) = 8^^ = ^ 

Note that according to (6.3), each function Kj^iG, e) has a discontinuity at g/ V2 but is 
analytic on each of the two intervals of the partition. Hence, it is piecewise analytic on 
(Q,g] in the sense of Section 8.5. 

From (6.1), the scale ^ appears with multiplicity 4*, for A; = 0, 1 , 2, . . . , so the scaling 



zeta function is 



1 



s e 



(6.5) 
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I I 



JULI 




*(UG,) *^(UG,) TtO) 

Figure 6.8. The U-shaped example U discussed in Section 6.2. 



T 

ee/i V ee/2 



T 




Figure 6.9. The generator G - G\ of U from Figure 6.8; see Section 6.2. 



It follows that the scaling complex dimensions are simple, and given by 

Dj:-{D + inp ■•. n € Z) with D - logj 4, p 
and the corresponding residues are 



In 
log 3' 



res (^(i); i = Z) + i«p^ 



1 



logs' 



for all n e ! 



(6.6) 



(6.7) 



Finally, we have the disjoint union - U {0, 1,2). All that remains is the substi- 
tution of the above quantities into the formula given in Corollary 5.9. We obtain 



(2 - k)K,iG) 



2-D-i«p 



z 

k=0 



Kk(G) 

1 - 4 ■ 3-* 



Jie) 



,2-t 



(6.8) 



where 7(e) :- max{j > 1 ■. dj^s < g) V as in (3.10), and [x] is the integer part of x. The 
computations for higher-dimensional analogues (like the Menger sponge) are extremely 
similar In each case, the only complication is to obtain the tube formula for the generator. 
Observe that 7~ is a lattice tiling in the sense of Proposition 5.5. 

6.2. U-shaped modification of the Sierpinski carpet. The U-shaped fractal of Figure 6.8 
is a modification of the Sierpinski carpet obtained by removing one contraction mapping 
from the self-similar system, and composing some of the remaining mappings with ro- 
tations of +7r/2. The generator G - Gi of U from Figure 6.8 provides an example of 
why it is useful to remove the requirement that lim£_40+ Kk{G, s) exists from Definition 2.1 
(Steiner-Uke)^; see Figure 6.9. 



^This assumption was part of the definition of Steiner-like in [LaPe3] but was removed in the present paper. 
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lid- — 



edge^, 



-E2 
Ei 



edge,, — . 

sector 
triangle 
solid 



Figure 6.10. The relation of V{G, s) to V{G, f ). 



To discuss G = Gi, let us consider the countable partition of [0,^) defined by the se- 
quence of intervals = 3^), for m = 1,2, Then the function m(s) :- [- log^ 2s] 

gives the index m - m(e) for which s e /,„. 

In this example, V{G,e) satisfies a recurrence relation (see the left-hand side of Fig- 
ure 6.10) given for e e /; by 

y(G,fi) = 9y(G,f)+^-^+j;r-^je2 for | < e < (6.9) 

The generator G is a union of countably many rectangles whose interiors are disjoint; 
consider these rectangles as defining a sequence of "chambers" {E,„}'^^^-^, as depicted in 
Figure 6.10. For s € /,„, the constant term in V(G, s) (corresponding to the region labelled 
"solid" in Figure 6.10) is given by the volume of U?lm+i which is 
00 \ 



\JE, 



\(=m+] 



1 „\ 1 11 



A2\—G\ = —A2iG)^ (6.10) 



6.3. A binary tree. In this section, we consider the example of a binary tree embedded 
in in a certain way. This example shows how a very slight modification can change a 
monophase generator to a pluriphase generator, and also how one can compute the tube for- 
mula for a set which is not a self-similar fractal (but which does have some self-similarity 
properties). 

Consider the fractal sprays depicted in Figure 6.1 1. Each of these figures is formed by 
an equilateral triangle whose top vertex is the point ^ - (1/2, V3/2) and whose base is 
the unit interval. Beginning at ^ and proceeding down one side of the triangle, one reaches 
the first branching at the point located | of the way to the bottom in (a) and at the point 
located | of the way to the bottom in (b). Consequently, the leaves of the first tree are the 
points of the usual ternary Cantor set, and the leaves of the second tree are the points of the 
(self-similar) Cantor set which is the attractor of the system {^P^-x) - |, *l'2(x) = | + |). It 
is clear from the "phase diagram" to the right of each spray that (a) is monophase and (b) 
is pluriphase. 

6.4. Apollonian packings. We consider the fractal spray associated to an Apollonian 
packing; see Figure 6.12. Recall that the construction of this packing begins with three 
mutually tangent circles contained in a disk which is mutually tangent to all three. For the 
next stage of the construction, a new circle is inserted into each lune so as to be tangent 
to its three neighbours. The Apollonian packing is obtained by iterating infinitely many 
times. After removing the outermost disk, the rest of the circles in the packing form a 
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(a) pj=3-' (b) pj=4-J 



Figure 6.11. The two binary trees discussed in Section 6.3. The left- 
hand one has a monophase generator and scaling ratios of the form £j = 
3 the right-hand figure has a pluriphase generator and scaling ratios of 
the form (j - A'K 




Figure 6. 12. The first three stages of the construction of the Apollonian 
packing for three circles with equal radii. The associated fractal spray 
does not include the outermost circle. 




Figure 6.13. The action of the Apollonian group on a configuration of 4 circles. 



fractal spray whose (monophase) generator is a disk, by [GLM^2, Theorem 4.1]. This ex- 
ample of a fractal spray was suggested to us by Hafedh Herichi. Full details on Apollonian 
packings and the Apollonian group may be found in [GLM^l, GLM^2]; we recommend 
the lecture notes [Sar] for an introduction. 

ApoUonius' Theorem states that given any three mutually tangent circles Ci,C2,C3, 
there are exactly two circles C4 , that are tangent to the other three (allowing the possi- 
bility of a straight line as a circle of infinite radius). Thus, if we have any configuration of 
four circles, one may be removed and replaced by its counterpart; see Figure 6.13. 
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Let a - (ai,a2, 03,04) be the 4-tuple whose entries are the reciprocal radii (i.e., the 
curvatures) of the four circles in a mutually tangent configuration. Descartes' Theorem 
states that these numbers must satisfy F{o) = 0, where F is the quadratic form 

F(o) — 2(flj + O2 + Oj + 04) - (fli + 02 + 02 + 04) . (6.11) 
If we start with circles of given radii Oj^^flj', 03 then this allows us to find a fourth via 

O4 — 0[ + O2 + Oi, + 2 y]Oi02 + 0\Ot, + 020^. 

Thus, if we start with three circles of radius Aj^' = a:^' = Oj' = 1, as in Figure 6.12, then 
the mutually tangent circle which encloses them will have radius ' = (3 - 2 V3)"' . 

For a starting configuration of four mutually tangent circles where one has negative cur- 
vature (so it encloses the other four), as in the top of Figure 6. 13, one can use \heApollonion 
group (a subgroup of SLuiT) generated by matrices A], A2, A3, A4) to geometrically obtain 
the other circles of the packing. Beginning with the configuration o - (01,02,03,04), one 
replaces a circle C,- with C- (its reflection with respect to the other three) and the new inra- 
dius is obtained from the corresponding matrix multiplication. For example, swapping the 
first circle Ci with its reflection C[ yields o' = (o'^, 02,0^,04) = {01,02, 03,04)Ai, where 
1 /a[ is the inradius of the new circle C'j . 

Consequently, the scaling zeta function ^ may be determined by collecting (with the 
proper multiplicities) the reciprocals of the entries of the 4-tuples 

00 

U U 

«=0 aieW" 

where W" = {1,2, 3,4)" is the collection of 4-ary words of length n, and A,y := A,,A„ ...A,_ 
is the matrix product corresponding to w = /1/2 . . . /„ e W". If G is a disk of radius r, then 
it has inner tube formula 

Vcis) ^ 27TrE - 7TE^ . (6.12) 

This example has a monophase generator with coefficients kq = -n and ki = 2n, when 
r = g=p(G) = 1. 

Given the radii of four mutually tangent circles (one of which contains the other three), 
it is possible (but nontrivial) to determine the radii of the other circles in the packing (and 
also to obtain asymptotics for their rate of decay; see [GLM+2, GLM^l]). Therefore, we 
omit the full details of the tube formula for Apollonian packings. Note, however, that up to 
a normalizing factor, the scaling zeta function of the Apollonian packing is given by ^(i) = 
YjJ=\ where X. = {^jljli denotes the sequence of radii (i.e., the reciprocal curvatures) 
of the circles comprising the packing (written in nonincreasing order and according to 
multiplicity). 

From [GLM^2, Theorem 4.2], the collection of circles in any packing is known to have 
Hausdorff dimension 

1.300197 <dimH < 1.314534. 

Note that all packings have the same Hausdorff dimension, since any two packings are 
equivalent by a Mobius transformation. Note also that the Hausdorff and Minkowski di- 
mensions coincide for Apollonian packings; see [Boy] (and compare [LvF2, Theorem 1.10]) 
We are hopeful that the methods of this paper will assist in the study of the Minkowski di- 
mension of such objects. Moreover, the determination of the complex dimensions of an 
Apollonian packing is an interesting and challenging problem. 
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7. Proofs of the main results 
7.1. Proof of the tube formula for fractal sprays, Theorem 4.1. 

The proof of Theorem 4. 1 is divided into several steps. We begin with a discussion intended 
to motivate and explain the approach. 

Remark 7.1 (The philosophy behind the head and the tail). A technical part of the proof 
of Theorem 4. 1 is inspired by the proof of the pointwise explicit formulas given in [LvF2, 
Chapter 5 and Section 8.1.1]. The main idea underlying the proof, however, is new and 
relies on the notions of "head" and "tail". 

In Section 7.1.1, we will split the tubular zeta function into a finite sum, which we call 
the head, and an infinite sum, which we call the tail (denoted by Wuiv respec- 

tively); see (7.4). Similarly, we will split the tube formula into corresponding finite and 
infinite sums (Vhead and Vtaii, respectively); see (7.15). This decomposition will allow us 
to avoid repeating the same type of argument as appears in the proof of Proposition 3.6, 
in several different instances. In particular, the decomposition into "head" and "tail" pro- 
vides a technical device which allows us to use the Heaviside function as expressed in 
Lemma 7.4. 

The "head" and "tail" decomposition is justified by the observation that for every fixed 
7 e N, the complex dimensions of a fractal string X - {^>)°li do not depend on the first 
J scaling ratios i^,. . This is the idea underpinning Proposition 3.6: the zeta function 
of the string X.j ■= X.\{(\,(2, ■ ■ ■ - K;)}ly+i has the same poles as the zeta function 
^ of the full string X.- Indeed, the function / defined by f{s) - p", where p > is some 
positive real number, is entire, and so is any finite sum of such functions. Hence, if ^ is 
meromorphic in some connected open neighborhood Q. of W, then since 

J 

for each s e Q, the truncated zeta function is meromorphic in Q and a point o) e D. 
is a pole of ^ if and only if it is a pole of . Moreover, the residues of ^ and at 
a> obviously coincide. The same applies if the first scaling ratios are changed instead of 
omitted. 

In addition to being technically useful, the "head and tail" decomposition helps one 
understand the conceptual difference between the contributions made to the tube formula 
by the integer and scaling dimensions, and the origin of the error term. Indeed, the bulk of 
the proof of Theorem 4.1 lies in showing that Vhead is given by the residues of ^fi^^.^iis, s) 
at the integer dimensions, and Vtaii is given by the residues of ^7-,aii(e, s) at the scaling 
dimensions; cf. (7.16) and (7.17). 

7.1.1. Splitting the tubular zeta function. For A: = 0, 1, . . . we define the function fk : 
(0, 00) ^ R by 

fk(e):^Kk(G,e)-K,(G). (7.1) 
This function measures the error of replacing the A:* coefficient function Kit(G, ■ ) by the 
constant Kk(G). Note that fkis) = for e > g and that in case Kk{G, ■ ) is constant, fk = 0. 

By employing the functions f^ defined above, the tubular zeta function can be rewrit- 
ten in a more convenient way (see Definition 3.5): 

j=l A=0 j=\ \k=0 J 
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In the second term, it is now possible to separate the scaling zeta function since the 
sum over k does not depend on j any longer. Moreover, according to the fact that (j — > 
as j — » oo, the first sum is taken only over finitely many integers j, namely, those indices 
for which (jg > e. (In fact, it is the finiteness of the first sum which ensures that the 
expression on the right-hand side in (7.2) converges absolutely exactly when the series 
defining ^j- does. Hence, (7.2) (as well as (7.4) below) holds for all e > and all 5 e C 
such that the second sum converges.) Set 

pj:^p{&)^tjg, (7.3) 

and recall that J - J{e) is the largest index such that pj > s, cf. (3.10). Obviously, pj is 
the inradius of the set - ^/G) and G^ c (G^)_£ iff ; > J. Hence, for all j > J, 

f]^e>8 and f,,({j^s) = for k = 0,1, . . . ,d. 

Interchanging the order of summation in the first term of (7.2), and making use of (2.7), 
we conclude that the tubular zeta function is given by 



k=o j=l 



u=o 



(7.4) 



with J - J(s), as in (3.10). Note that the t/* term in ^ftmjis, s) vanishes because of the 
presence of the factor (d - k) inside the sum. 

In combination with the splitting in (7.4), Proposition 3.6 yields the following result. 

Theorem 7.2. Assume that W is a window for ^ and that Q is a connected open neigh- 
borhood ofW in which ^ is meromorphic. Fix an arbitrary e > 0. Then: 

1 (Meromorphic continuation and poles of ^j-). Both i^j-{e, ■ ) and ^7-,aii(e, ■ ) are 
meromorphic in all of Q. Furthermore, the set of visible poles of these functions is con- 
tained in DfiW) — Dt ^> given in Definition 3.7. 

2 (Head and tail decomposition of i^j-). For every s e Q. (in particular, for every 
s e W), tlie meromorphic continuation of Cf^ie, ■) to D. is given by 

^^(s, s) = ^r,h.J^^ + ^T;.aii(e' ^) • (7.5) 
Here, i^r.h,J>^^ ' ) given by 

k=o " M 

which is meromorphic in all of C, with poles in {0, \,. . . ,d], where fk is as in (7.1) and 
fjig ^ „i X/ n5 

^T;taii(e' ■)toQ. is given by 



7(e) := maxjy > 1 ■. {j s < g)y 0, as in (3.10), while the meromorphic continuation of 



/:=() 

where i^£(s) denotes the meromorphic continuation ofYIJ^i {"j, os usual. 



^It may be useful to keep in mind that even tliough 7 = 7(e) is finite for every e > 0, it tends monotonically 
to oo as £ — > 0^, since tj decreases monotonically to 0. 
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3 (Residues of ^^). For each u € DriW) = (©x U {0, 1 , . . . , t/)) n W, we have 

res (^^(e, s); i = w) = res (^7;h,,d(e, s); s^u)+ res (^7;„ii(e, s); s = , (7.8) 
which implies 

res (^^(e, i); 5 = = res (^7^,„,(fi, s); s = (7.9) 

for each u e \ (0, 1, ■ • ■ , ^Z}, iwce ^7-head(^' holomorphic on C\ {0,1, ... ,d}. 

(i) In the case when to G D^iW) \ {0,1, . . . ,d} is a simple pole of ^^{s, ■ ) (and hence also 
a simple pole of^j-), 

res {^^{e, s); u) = res (^^(i); u) V ^(d - k)Kk(G). (7.10) 

A:=0 

f»'j In the case when o) - k & {0,1, . . . ,d], we have 

m 

res {(r,,.Je^ k) = e'-' J] ■ (7-11) 

7=1 

Furthermore, if cd - k & {0,1, ... ,d — \] is a simple pole of^j-ie, ■ ), we also have 

res {^^^Js, s); k) = s'-'^^ikWC). (7.12) 

If a> — d € W and D < d (as assumed in Theorem 4.1), then d is not a pole of ^j.. Hence, 
formulas (7.7) and (2.7) imply that 

res (^^,^„(£, s); d) = ^^{d) (Kd{G) - A,{G)) . (7.13) 

We leave it to the interested reader to perform the necessary (and elementary) com- 
putations needed to deal with the case when w is a multiple pole of see [LvF2, Sec- 
tion 6.1.1]. 

Remark 7.3. If we assume that oj e {0,1, . . . ,d - 1) n DjiW) in part 3 of Theorem 7.2, 
then one can express res (^7;,^ii(e, s);s - oj^ in terms of the constant term of the Laurent 
expansion of ^(s) at s = A; (even in the case when is a multiple pole of ^), much as was 
done when d = I 'm [LvF2, Corollary 8.10 and Remark 8.11]. 

7.1.2. Splitting the parallel volume. For fixed e > 0, we split the inner parallel volume 
V{T, s) in a similar way as the tubular zeta function above. Taking into account first (2.5) 
and then (7.1), we can use (7.3) to write 

j:pj>e j:pj<B 
d 

= Z Y,^''-''(]Kk(G,rj's)+ ^ (7.14) 

i:Pi>E k=0 r.pj<£ 

j:pj>B k=0 j-Pj>E k=0 J-Pj^'^ 

Recall that the sum over j in the first two terms is finite for each fixed e > and that 
the number of terms is given by 7 = /(s). Therefore, in both terms, the sums can be 
interchanged. In the third term, the homogeneity of the (i-dimensional Lebesgue measure 
implies Ad(G^) - AditjG) - C'jAdiG) for each j. Thus (2.7) yields 
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k=0 j=l k=0 j-Pj>E fPj^E 

k=0 7=1 k=0 \fPj>E j'Pji" 

=:Vhead('r,£) =:V,„ii(r,e) 



(7.15) 



7.1.3. Outline of the remainder of the proof of the tube formula. In light of (7.15), the tube 
formula (4.1) of Theorem 4.1 will follow upon verification of the following two assertions: 

Vhead(T, e) = ^ res [^r,h<^Jji, s); i = + '??head(e) , and (7.16) 
Vtaii(T-, e) = ^ res (^^,^,(e, *); s^cJ) + A,{G)i^j.{d) + KM, (7.17) 
where the error terms "Rhead and ^taii are given by 

[5(0)1 

^head(e) = ^ res [CT,h„ Ji^^ ^)'^ s ^ k) , and (7.18) 

il:=0 

^taii(£) = ^ ^ 5r;,a,i(«' ^)ds , (7.19) 

respectively. (Here [x] denotes the integer part of x.) Indeed, the assumption 5 (0) < 
implies immediately "Rhead = and, by (4.2), we have 'R{s) = 9?taii(e)- Therefore, if the 
formulas (7.16) and (7.17) hold, then since ^7-(e, s) = {j-^i,^^iie, s) + ^j-^^^iie, s) by (7.4), it 
follows from (7.15) that 

y(r,e) = yhead(r,e) + y,„i(r,e) 

= ^ res [^r^Js, s); s ^ oj) + Keadie) 

+ res [^r,vJs, s); s = w) + Ad{G)(j:{d) + 
= Yj ""^s {^rM..M ^) + ^T;uui(e' s = + AAG)^j-(d) + Hie) 
= Y res '^); * = + ^rf(G)^x(^) + ^(«) ' 

and we obtain the tube formula (4.1), with IHe) given as in formula (4.2): 

Kis) = KM = ^ ^ (7-20) 

Observe that in light of (7.4), the integrand Crt^J^^^ coincides with that of (4.2). 

We note that in order to complete the proof of Theorem 4. 1 , we will still have to establish 
that the error term '??taii(e) (which coincides with the full error term Kie), since '??head(fi) = 0, 
as noted in (7.20)) satisfies the asymptotic estimate '??taii(e) = (9(e''"™P'^) as e — > 0^ (in the 
languid case) and that "RtaiKe) - for all < e < vnm{g,A^^g] (in the strongly languid 
case). This will be accomplished, respectively, in Section 7.1.7 and Section 7.1.8. 
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We will establish independently (in Section 7.1.4 and Section 7.1.6, respectively) that 
VheadC^.fi) can be expressed as a sum of residues of ^7-i,(.„j(e, s) (as given by (7.16)) and 
Viai\{T, s) as a sum of residues of (ftmii^' (^^ given by (7.17)). While ^7-|,ead(^' ' ) mero- 
morphic in C and has poles only at the integer values 0,1, ... ,d, the function ^7-,ij||(e, s) 
recovers the full set of the complex dimensions of X. (compare Theorem 7.2). This fits per- 
fectly with the observation (compare Remark 7.1) that the first lengths of a fractal string 
(or, as here, the first scaled copies of G in a fractal spray) do not affect its complex dimen- 
sions. The scaling complex dimensions are hidden in the "tail". The derivation of the first 
part is elementary and exact in that the error term just collects the residues at the integer 
dimensions not contained in W. The derivation of the second part is more involved. It uses 
techniques similar to those used in the proof of [LvF2, Theorem 8.7], the pointwise tube 
formula for ( 1 -dimensional) fractal strings. Here, it is necessary for X to be languid (for 
the first part of Theorem 4. 1) or strongly languid (for the second part). 

7.1.4. Proof of (7.16). The residue of ^7-head(*^' ^'^ - ^ ^ {0, 1, ...,£/) is given by 
formula (7.11) in Theorem 7.2. Observing that the k-th term of VheadC^, e) in (7.15) has 
exactly the same expression, we conclude that 

d 

Vhead(T, e) = ^ res [^T^t^^Js, s); s = k) . 

k=0 

For a fixed screen S (and a corresponding window W) such that S(0) i {0, . . . , d], we can 
split this sum into two parts, according to whether k is contained in the interior of W or in 
the complement W . Since (j-,^^.J(E, ■ ) has no poles outside the set {0, . . . , d], we can safely 
extend the first sum to include the residues at all complex dimensions visible in W. Thus 

Vhead(T, e) = ^ res [^r,h.;M^ s);s^a))+ res [^r,h. J>^^ s ^ k) , 

where the second sum is !Rhead(e) (as defined in (7.18)). This completes the proof of (7.16). 
We note that in the first part of Theorem 4.1, the assumption 5(0) < for the screen 
ensures that "Rhead = 0. In Corollary 4.2, where this assumption is dropped, we also have 
^head = 0, but for a different reason; indeed, in the monophase case, i^fhc'^ii^^ s) itself 
vanishes identically. In the general case (which is not treated in this paper), when the 
integer dimensions may not all be visible (i.e., 5(0) > 0) and when G is not necessarily 
monophase, the error term '??head will have to be taken into account; see also Section 8.4. 

7.1.5. The Heaviside function. Before continuing on to the proof of (7.17), we need to 
make some remarks on a certain useful form of the Heaviside function // : R — > {0, 1), 
which is defined (as in [LvF2, (5.10)]) by 



0, jc < 

1, x = 0, (7.21) 
1, X > 0. 



Following the proof of [LvF2, Theorem 8.7], we will exploit the following integral 
representation of the Heaviside function, which comes from number theory [Dav, p. 105], 
and was refined in [LvF2, Lemma 5.1]. 

Lemma 7.4. For x,y,c > 0, the Heaviside function is given by 

1 f^'" ds 
H(x-y)^— xY'-. (7.22) 

27ri Jc-ico s 
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H(x-y) 



s 



J- > 0, the 
(7.23) 



1 Proof. Taking A; = 1 in [LvF2, Lemma 5.1] for a screen with T - T+ 

2 Heaviside function is approximated by 

where the absolute value of the error term E is bounded by 

{{xy'^)' Y min{r, I log(;cy-i)r 1 ), ifx+y. 
If, \ix-y. 
It is now easily seen, that, for arbitrary fixed values x,y,c > 0, the error term E - E{T) 
vanishes as T — > oo. Since this is true for all x,y,c > 0, the result follows. □ 



(7.24) 



6 7.1.6. Proof of (7.17). The proof of this part follows roughly the lines of the proof of 

7 [LvF2, Theorem 8.7]. One can rewrite the expression of VtmiiT, s) in (7. 15) as 

d-i f 

Vtaii(r,e) = 2g-Vi(G) Yj X p1 + i^diGX^id). (7.25) 

8 For k - 0,1, . . . ,d - 1, denote the expression within the parentheses of (7.25) by Vi(e). 

9 Using the Heaviside function as defined in (7.21), we write 



vde) = 2 ""'''P' ^YuP'j^Y. [«'"VJ^(P; pp(s - pj)] . (7.26) 

j-Pj^e j:pj<E j=l 

10 Note that, in case s = pj for some j, the corresponding J* term in the sum on the right-hand 

11 side of (7.26) is 

e'^VjHipj -s)+ pp(s - pj) = ipf + Ip'j = P'p 

which equals the value given by the left-hand side of (7.26). 

Now, fix some constant c such that d-l < c < d and D < c. (This is possible, since 
the abscissa of convergence D of ^ was assumed to be strictly less than d.) Then c-k and 
d - c are positive numbers and so, by Lemma 7.4, the j-th term in the above sum is given 
by 

r^^'^s'p".-'^. 

Substituting s = t + kin the first integral and s = d - tin the second one, we get 

. Pj'' 



12 
13 
14 
15 
16 



■C-k + ioO J, 

d-k-tk+t 



pc-k+i 
*Jc-k-ic< 

'( in the 1 

. ^''PJ 

«7C — lOO 



' 2OT(i - k) Jc-ica 2jTi(d - s) 

Combining these integrals, one sees that (7.26) can be rewritten as follows: 



fc p. 



*^ ^ 27T j-^^ (c-k + it){d - c - if) 



dt. 



(7.27) 



19 
20 



This last integral converges absolutely, which is seen as follows. Since c was chosen 
strictly between d-l and d, the numbers d - c and c-k are positive and so we have for 



21 all real numbers f. 



\(c-k + it)(d - c- it)\ = |(c - k){d - c) + + i(d - k)t\ > (c - k)(d - c) + f-. 



22 Hence 



t p. 



(c - k + it){d - c - if) 



8^-'p)\s-"\\p';\ 



e'^-'p'. 



\(c-k + it){d -c - if)| (c - k)(d - c) + ' 
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which impHes 



Ly f — f — P_i dt<^yp^/-^ f 



dt 



.^^ ^-^{c-k){d-c) + fi 

1 „ , , r" dt 



1 r' 

^ Jo 



(c - k){d - c) + f 

Note that the last expression is finite because c was chosen such that c > D, where D 
is the abscissa of convergence of It follows that the integrand in (7.27) is absolutely 
integrable and it is safe to interchange the order of summation and integration. Hence 

where the sum YjJ=i P] can be replaced by g''^j.{s). Inserting the derived expressions for 
Vf:{s) into (7.25), we obtain 

1 ' nc+ica 2''^^(d — k) 

Vtaii(T-, e)^ — y K,(G) s''-'i^(s) /^ ; ' ds + Kd{G)C^{d) 

—7^^y^{d-k)Kt(G)ds + KAGX^id) 

= z-^ ( ^rtaii(*^' '^■s + Kd{G)^j~{d), 

Zm Jc-ioo 

by (7.4). Since li is not a pole of it is a simple pole of ^7-,;,ii(e, ■ ) and (7.13) yields 

res {^^^e, s); s ^ d) ^ ^j:(d){KAG) - Ad{G)). (7.28) 

Therefore, we obtain 

Vtaii(T, ■s) = ^ J (tmi(^^ *)^'' + res (^T;ia,i(e' '^); = t/) + ^_c(d)Ad(G). 

Now the machinery of the Residue Theorem can be applied. When pushing the line 
of integration towards the screen S, we collect on the way the residues of the poles of 
^7~iaii(^' ■ ) '^hat lie between the line Re i = c and S (see the proof of [LvF2, Theorem 8.7]). 
The definition of the screen S and the window W imply that ^ is meromorphic in W 
and, since D < c, ^j. has no poles to the right of the vertical line Re i = c. Therefore, 
by Theorem 7.2, any pole of ^7-,i,i|(e, •) in the region between Re i = c and S is either 
contained in {0,1, . . . ,d - 1) or a pole of ^ in W, i.e., an element of DTiW) \ {d}. Recall 
that ^7-,aii(e, ■ ) has another pole at d but, since d is not passed when pushing the line of 
integration towards the screen, it does not occur again. 

At this point, the languidness of ^ comes into play. Using the sequence {T„]nez of 
Definition 3.3, we write VtaiiC^, e) as a limit of truncated integrals: 

I rc+iT„ 

VUT, s) = Hm — ^^^Js, s)ds + res U^^Js, s); s ^ d) + ^^(d)AAG). (7.29) 

If we now replace the vertical line segment C\„ := [c + ir_„,c + iT„] of integration by the 
curve given by the union of the two horizontal line segments and the truncated screen S\n, 
that is. 



U\„ := [c + iT„,S(T„) + iT„], 
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(\„ := [c + ir_,„5(r_„) + ir_„], and 
5|„:={5(f) + if:fe[r-,„r„]) 

with proper orientations, the Residue Theorem impUes that the n* integral in (7.29) is 
equal to 

X {ir,vJ'^^ s)\s^Lj)+ %,{s) + ul (e) + lI (s), 

where D(W\„) := Dt-(W) n W\„ is the set of (possible) poles of ^7-,^ii(e, ■ ) that lie inside 
the region W\„ bounded by the curves U\„,S\„,i\„ and C\„. Hence, W\„ is the "truncated 
window" associated to the truncated screen S\„. The term '??„(e) is given by the integral 

%As) = ^ ^T;,a,i(fi' ^) ds (7.30) 

and U„ (e) and £„ (s) are the corresponding integrals over the segments U\„ and £\„, respec- 
tively (traversed clockwise around W\„). More precisely, uI,(e) is given by 

(e) = z-^ ^r;taii(s, i) 

Js(7'„)+ir„ 



= J,(,„) + ^^"^ 2. (f _ ^ + iT„)(d - 1 - iT„r 

and is absolutely bounded as follows: 

If./. r , -t^mv ig'-v(G)i(af-^) 

According to the languidness condition LI of Definition 3.3 and the hypotheses of the first 
part of Theorem 4.1, there exist constants C > and y < 1, such that |^(? + ir„)| < C{T„ + 
1)''. Moreover, \t - k + ir„| > T„, for all A: = 0, . . . , li - 1 and, similarly, \d - t - ir„| > r„. 
Hence we get 

|t/l(.)l<:^C(r„ + l)^y^^^Mglz^ r (7.31) 

Since S{T„) > inf 5, the integral in this expression is bounded by a constant independent 
of n. Thus, there is a constant C\ > 0, independent of n, such that 

l^i(g)l <Ci " 2 ■ (7.32) 

n 

With similar arguments, one can show that the integral („ (e) is absolutely bounded by 
C2\T-nY'^(\T-n\ + \y, for some constant C2 > independent of n. If we now take limits as 
« — > 00, then T„ 00 and r_„ — » -00. Since y < 1, this implies \ ul(e)\ and (e)| tend to 
as n — > 00. 

7.1.7. Estimating the error term. To complete the proof of formula (7.17), it remains to 
show that the limit '??taii(£) '■- linin^oo '^n(e) exists and satisfies the asymptotic estimate 
■^taiKe) = (9(e''"*'"P'^) as e — > 0^, for which we utilize assumption L2 of Definition 3.3. 
Recall from (7.30) and (7.4) that the integral %j{e) is given by 



2ot Js|„ ^ {s-k)(d-s) 
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= — e Cc('^ W + lO > , 7777-— t; — 777:1 — 777-; — —{S it) + i)dt, 

where 5 '(f) denotes the derivative of S at t. Note that, since S was assumed in Defini- 
tion 3.2 to be Lipschitz continuous with constant Lip^, 5"(f) exists for almost all f e R 
and \S'{t)\ < LipS at those points. Hence the integral above is well defined and absolutely 
integrable, which is seen as follows: 



— E %£iS {t) + it)\ ) -— — . — \S (t) + i\dt 

M(e)(l + Lip5) ^ r^" ICr(5(0 + iOI 



y \Kk{G)\{d - k) \ ^ , (7.33) 

2n U ' Jt_„ \S(t)-k + it\\d-Sit)+it\' 

where the number M(e), defined by 

Mis) = max{e''-™P^e''-'"") ■ max{g^"P^g"") ■ max{l,g-''), 

is a uniform upper bound (in f) for the term s'^^^^'^g^^'^^'', for k - Q, . . . ,d. Now we use the 
languidness assumption L2, which states that there exist constants C > and y < 1 such 
that 1^^(5(0 + ir)l < C|r|^ for all |f| > 1. Observe that, since the screen S avoids the poles 
of the expression |fs(5(f) + if)l is bounded on any finite interval for t. Therefore, L2 
is equivalent to assuming that there are Ci > and 7 < 1 such that 1^5(5 (f) + if)l < Ci|fP 
for all \t\ > fo, where to is some arbitrary but fixed positive constant. (Simply choose C\ 
sufficiently large.) Next, we describe how to choose fo- Since the screen S is assumed to be 
Lipschitz continuous and to avoid the numbers {0,. . .,d} when passing the real axis, one 
can find positive constants fo and ro such that \k- S (f)| > ro for all |f| < fo and k - 0, . . . ,d. 
(That is, in a tube of width fo around the real axis, the screen S has at least distance ro to 
any of the lines Res - k, for k = 0, . . . ,d.) 

Now, for the remaining integrals in the above expression (and n sufficiently large), we 
split the interval of integration (r_„, T„) into (r_„, -fo) U (-fo, fo) U (fo, T„). In the first and 
the third intervals, we use (the modified) condition L2 and, furthermore, that 1^-5 (f)-i-if| > 
|f| and \S(t) - A; + if| > |f| to see that, for = 0, 1, . . . , - 1, 

jT_„ \s{f)-k + it\\d-s{t) + it\ - 'Jt_„ ° ' 

and, similarly, that the ^* integral over the interval (fo, T„) is bounded by the constant 
r-i V " /■ 

In the interval (-fo, fo), |{^(5(f) + if)| is bounded by a constant, say M, \S(t) - A: + if| > 
\S(t)-k\ > ro and, similarly, \d-S(t)+it\ > \d-S(t)\ > kq. Therefore, for^ = 0, 1, . . . 



X'o 
1(1 



|^^(5(f) + if)| 2Mfo ^ 

dt < — r— =: C2. 



\S(t)-k + it\\d-S(t) + it\ ~ rl 



Observe that the derived estimates for the A:* integrals are independent of k. Thus, putting 
the pieces back together, we have that (7.33) is bounded above by 

- tV - \T-nV-') + Czj , (7.34) 

where C(e) :- '^^^''>'-\^^^'('^'> YflZo \'^k{G)\(d-k). Consequently, ^„(e) is absolutely integrable 
for each n and e > 0. Moreover, since (7.34) converges to some finite value as « — » 00 
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(because y < 1), it follows that also ^taii(e) is absolutely integrable and thus integrable; 
i.e., '??taii(£) is finite for each £ > 0. Hence, the error term "RtaiKe) is given as claimed in 
(7.19). Finally, note that W|„ ^ W n {Re 5 < c) and 'DriW) n {Re > c) = [d] imply 
£)(W|„) Dt(W) \ {d}. This completes the proof of formula (7.17). 

Furthermore, from (7.34) and the definition of M{s) (see the discussion following (7.33)), 
it is clear that there is a constant C > such that |'Rtaii(e)l < Ce^'^"'''^ for all < e < ^; 
i.e., Hiaiiis) is of order 0(£^'^"p^) as e — > 0^. Recalling that H - Hf^i, this completes the 
proof of the languid case in Theorem 4. 1 . 

7.1.8. The strongly languid case. Now assume that ^ is strongly languid of order y < 2, 
as in Definition 3.4 and the second part of Theorem 4.1. Then there exists a sequence S„, 
of screens and corresponding windows W„, with sup5„, — > -oo as m — > oo such that LI 
and L2' are satisfied for each m (with constants C,A > independent of m). In addition, 
we may assume without loss of generality that sup 5,,, < S{Q) for all m > 1 (see the 
discussion preceding Corollary 4.2). For each screen S„, and for fixed « € N, consider the 
truncated screen 5m|„ (truncated at r_„ and T„) and the corresponding truncated window 
W,„\„ bounded from above and below by the horizontal lines Im s = T„ and Im s = r_„ and 
from the right by the line Re i = c. By the Residue Theorem, for each m and «, the n* 
integral in the sequence of truncated integrals in (the counterpart of) (7.29) is given by 

X 1"^^ (^7;.aii(e' i = w) + %n\n(s) + f/,{|„(e) + L{|„(e), 

ajeV(W„,,„) 

! f 

iust as in the languid case, and the integrals U , (e) and L , (b) over the horizontal line 
segments are similar to U„ (s) and L„ (s) above, with S replaced by S,„. First we keep n 
fixed and show that 'Rm\n{s) vanishes as m — > oo. Note that 'Rm\n{s) is given by the same 
expression as "R„(e) in (7.30), except that the integral is now over Sm\„ instead of S\„. Its 
absolute value is bounded by 

- E "'^ '|^^(5,„(f) + if)l >, ttttt; — , ^ .,11 , — rTTTT^I'^'BW + 'Iflf 
27r Jt ^ j;^\S„,{t)-k + it\\d-S,„{t) + it\ 



A=0 



B+1 r^"i, ^ ,A tl4'r(5„,(f) + 

yk,(G)i(j-^) ^^-^..cy.o-^^^^y — 'Idt, 



2- ..0 

where we used the inequality \S mit)-k+it\\d-S ,„(t) + it\ > |fp. Moreover, we utilized that, 
since the functions S „, are assumed to be Lipschitz continuous with a uniform Lipschitz 
bound B = sup^jLipS',,, < oo, the inequality \S'^{t) + i\ < B + I holds, whenever ^^^(f) 
is defined (which is the case for almost all f e R, independently of m). Now, by L2' of 
Definition 3.4, there are constants A, C > 0, independent of n and m, such that, for all f e R 
and all m e N, \^j.(Si„(t) + if)| < CA^^"'^'^^{\t\ + 1)''. Therefore, there exists a constant Ci, 
independent of n and m, such that 

For m sufficiently large (indeed, without loss of generality, for all m > 1), we have 5,„(f) < 
and so -Sm(t) - |5'm(f)l- Thus, provided that e < A'^g, we can bound the expression 
(g/^)-5,„(0^|5„(0l = (sA/gf-^'^^ from above by (£A/g)l'"P^'"l, which is independent of t and 
can thus be taken out of the integral. The remaining integral has a finite value for each n. 
Letting now m — > oo, | sup 5^1 oo and so |^„,|,!(£)l vanishes. 
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f 

When taking the limit as m — > oo, the expression f/;„|„(e) extends to an integral over the 
whole half-line (-00 + iT„,c + iT„] and Lj^^Js) to an integral over (-00 + iT_„,c + iT_„]. 
More precisely, uj^is) lim,„^oo uj^^^is) is given by 

2711 J _^ ^ {-f(t-k + iT„)(d-t-iT„) 



k=0 

By exploiting the languidness condition LI (which now holds for all f e R) and the in- 
equalities \t-k + iT„\ > T„ (for k - 0, . . . ,d - 1) and \d-t-iT„\ > T„, it is easily seen that 

C 

there exists some constant C2 > 0, independent of n and m, such that t/|„(e) is absolutely 
bounded as follows: 



\K^^C2'-^^£[^-j dt. (7.35) 

The remaining integral is finite, provided that s < g. Now, as n — > 00, |f/|{(e)| vanishes, for 
each e < g, and with completely analogous arguments, the same can be shown for |f/||^(e)|. 
Hence the tail volume is given in the strongly languid case by 

i.e., equation (7.17) holds without error term for s < mm{g,A^^g}. This completes the 
proof for the strongly languid case, and thus of all of Theorem 4. 1 . 

7.1.9. Proof of Corollary 4.2. In the monophase case, we have fk(s) = for all e > and 
for each k - 0,1,. . .,d. Therefore, (given in (7.6)) vanishes identically, implying 

^?head = (by (7.18)) and thus yhead(T,e) = for each e > (by (7.16)). Consequently, 
by (7.15), V(T, ■ ) = VtaiiCT, ■ ). Since -Rtai^e) = 'R{s), cf. (4.2) and (7.19), the assertion 
of Corollary 4.2 follows by observing that the assumption 5(0) < is not used in the 
proof of (7.17) given in Section 7.1.6 and Section 7.1.7. It is only used to ensure that 
the screen S avoids the integer dimensions 0,1,. . . ,d. Note that ('j-hem = also implies 
^j- - Hence, the error term 'R{s) is equivalently given by the integral (4.3) in this 

case, as explained in Remark 4.3. □ 

7.2. Proof of the fractal tube formula, Corollary 5.9. 

Before proceeding, we need to compute some residues. To this end, we introduce the 
tubular zeta function for the generator G. In addition to being a useful technical device, it 
reveals the structure of the residues of the tubular zeta function 

Definition 7.5. Let ^c(£j denote the tubular zeta function of the generator G, where G 
is assumed to have a Steiner-like representation as in (2.2). It is defined exactly as ^j-{e, s), 
except that the associated fractal string is given by with £\ - 1 and Ij - for 

all j > 2. In other words, it is the tubular zeta function of the trivial fractal spray with 
generator G. 

Exactly as in (7.4), we write - ^c.hcad + ^c.taii, so that for s € C, we have 

= If (737, 

(0, s>g, 

with fkie) - Kk{G, e) - K/^iG) defined as in (7. 1) for A; = 0, 1 , . . . , c/, and 
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^c,aii(£, s) = M,(G), e > 0, (7.38) 

a - s 

where 

M,(G) := - (7-39) 

To see why the second case of (7.37) should be true, consider the definition of fc.hcad in the 
counterpart of (7.4) and suppose we define 

Jc(e)■■^X^o,g)(E), (7.40) 
which is in parallel to (3.10), upon inspection. 

Observe that for every s > 0, (ci^, ■ ) is meromorphic in all of C, with poles con- 
tained in {0, 1, ... , d}. Hence, the set of "complex dimensions" of G consists of the integer 
dimensions {0, 1, . . . , li), and all of these poles are simple. 

Lemma 7.6 (Residues of ^c)- For < e < g, we have the following residues of^c- 

res (^G,head(£, s); s ^ k) ^ s'-'Ms), k^O,l,...,d, (7.41) 

res (^G,,aii(e, s); s ^ k) ^ e''"*^/ft(G), k ^ 0,1, . . . ,d - I, (7.42) 

res (^G,,au(e, s); s ^ k) ^ k^G) - A,(G), k^d. (7.43) 

Proof. In light of (7.37)-(7.39), each of ^G,hcad(£, ■ ), ^G,taii(£, ■ ) and (c(.s, ■ ) is meromor- 
phic in all of C, with (simple) poles contained in {0, 1, . . . ,d], for Q < e < g. To show 
(7.43), simply use (7.38) and (7.39) to compute 

d-\ 

res (^G.iaii(e, s); s ^ d) ^ lim {s - d)^CMiis, «) = - V, g'^^'^'^kiG) = kj{G) - A,i{G), 

using (2.7) to reach the last equality. □ 

The following result will not be used in the sequel but may be helpful for the reader; it 
provides a "residue formulation" of the given Steiner-Uke representation of G. 

Corollary 7.7 (Exact tube formula for G). For all e e (0, g], 

d 

y(G, e) = ^ res (ic(E, s); s ^ k) + A,{G). (7.44) 

k=Q 

Proof. First, note that it follows from (7.41) and (7.1) that for Q < s < g, 

d d 

^ res (^G,head(e, s);s^k)^^ e"'"* (/c^(G, e) - Kk{G)) 

k=0 k=0 

d d 

^YjEd-'K,iG,E)-Y,e'-'K,(G) 

k=0 k=0 
d-\ 

= V(G, e) - ^ E'-'K.iG) - kAG), (7.45) 

k=i) 

where we have used (2.2) in the last equality. Furthermore, by (7.42) and (7.43), we have 
fore e iQ,g], 



38 



MICHEL L. LAPIDUS, ERIN P. I. PEARSE, AND STEFFEN WINTER 



2 res i^cM^is, sy,s^k)^Yj ^''"'^kiG) + MG) - Ad(G)) . (7.46) 

Since = ^Chead + ^ciaib the result now follows by adding (7.45) and (7.46). □ 

As an alternative proof of Corollary 7.7, one can obtain (7.44) by applying the second 
part of Theorem 4. 1 to the trivial fractal spray on G. However, we feel that the proof given 
above is more edifying and more straightforward. 

7.2.1. The residues of (j-. Let T he a self-similar tiling with a fractal string X. - {fj}J^i 
and a single generator G for which a Steiner-like representation has been fixed. Let ^j- - 
i^^ie, s) denote the tubular zeta function of T, and let (j- - CT,h^ii + '^tm\ head-tail 
decomposition, as in Section 7.1.1. In light of (7.4), we deduce from (7.38) and (7.39) that 
factors as follows: 

^r,.,{e,s)^^cU^,s)Cj:(s). (7.47) 

Furthermore, still by (7.4), 

^T,...M^ = Z 4t l! (7.48) 

with fk as in (7.1) and 7(e) as in (3.10), as usual. Recall that 7(e) — > oo monotonically as 
e — » O'*", since £j decreases monotonically to as j — > oo. 

Remark 7.8. Observe that (7.48) is not at all the counterpart of the factorization given in 
(7.47). Indeed, it clearly does not enable us to write ^7-^^^,^ as the product of ^c.head and ^ 
(which would be false). This is the source of some difficulty if we wish to estimate the 
residues of ^7-|,ead(^' *) as e — > 0^. 

Lemma 7.9 (Residues of ^^). Fix s e iO,G]. Then: 

(i) When oj € T)^ . . . ,d] is a simple pole of^^, the residue res (^7-(e, s);s - oS) 
is given by 

^G,taii(£> w)res {^jis); s-uj^- res [^j-{s)-s^uj)m^{G), (7.49) 

with MJG) = Zto - kWG) as in (7.39). 

(ii) For o) — k e {0,1, . . . ,d} \ ©x,^ the residue res (^7;taii(*^' = ^) is given by 

I, , , A U'-'K^iGX^ik), k + d, 

res <T-t,ii(e> s)\s-k]-< (7.50) 
\{Kd{G)-Aa{G))^j~(d), k^d. 

(iii) For u) — k & . . . ,d), the residue res [(.-j-^^^^Jis, s); s — k^ is given by 

m 

res (^7;h.ad(e. s); s ^ k) ^ s"-' ^ {]fk{{fe), (7.5 1) 

with J(e) as in (3.10) and fk as in (7.1). 



^Note that if 7" is a self-similar tiling, then the only real pole of ^ is D = D£ < d. Hence, the only way that 
0) could belong to both T)£ and (0, 1 , . . . , rf) would he if oj = D = k, for some ie{0, l,...,d-l). 
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Proof. In light of the factorization formula (7A7), (7.49) follows from (7.38) and the fact 
that, under the assumption of (i), 

res {^.j-is, s); s = w) = res (^^j-^t^ais, s); s - oj^ , for cj e D£ \ {0, 1, . . . , li), 

while (7.50) follows from (7.42)-(7.43) of Lemma 7.6. Note that w is a simple pole of ^ 
in part (i), and hence it is at most a simple pole of ^fti,n(^' ■^)' whence 

res {(r,un^s, s); i = w) = lim(i - w)^^,^;/^' = ^G,iaii(e, w) res [^j^is); i = w) , (7.52) 

from which (7.49) follows in light of (7.38). Since w = A: e {0, 1, . . . , t/) is a simple pole of 

res {^tm J^' «); = lim (s - kXr^ J^' «)' 

and hence (7.51) follows immediately from (7.48). Finally, as was already observed, the 
poles of and belong to {0,1, ... ,d} and D7-, respectively. □ 

Remark 7.10. Note that Lemma 7.9 is valid for an arbitrary fractal spray satisfying the 
hypotheses of the first part of Theorem 4.1, but without the assumption that 5(0) < 
(which is not necessary for Lemma 7.9 to hold). 

7.2.2. The proof of Corollary 5.9. Let T he a self-similar tiling satisfying the hypotheses 
of Corollary 5.9. Note that since the poles of are assumed to be simple, it follows that 
and {0,1, ... ,d} are disjoint; that is, all poles of ^ are simple and D i {I, . . . ,d - 1). 
Recall that since 7~ is a self-similar tiling, we have < D < d and D is the only pole of 
^ on the real axis. See footnote 6. The following proof makes use of the decomposition 
= ^T;hcad + ^rtaii from (7.4). Since = D£'J{0,1, . . . ,d] is a disjoint union, the present 
hypotheses and Theorem 5.7 yield (for s e iO,g) and with e,t(e) defined as in (5.21), for 
k^O,l,...,d) 

V(T, e)^ Yj res {^tmM s); u) + ^ res (^^,,„(e, s); k) 

we»x ke{0,l,...,d} 



+ 2 ves{^^^,^Js,sy,k) + AAGX^(d) (7.53) 

ke{0,l,...,d} 



+ Y ^'''^kis) + UG)^£{d), (7.54) 

k=0 

from which (5.18) follows. In (7.54), we have set 

d-l w-k 

MUG)^y^(d-k)Kk(G) 

as in (7.39), so that c„ - M^(G) res (^(i); s - coj /(d - o)) for co e D£, and applied 
Lemma 7.9 to obtain the precise values of the residues of ^^^^^^j and ^7-,^,!,; see (7.49)-(7.5 1). 
In particular, this verifies (5.19). Note also that the residue of ^7-,;,ii(e, s) at s = d and the 
term Ad{G)i^j^{d) have been combined to yield Kd(G)i^j-{d) = cj. This verifies (5.20). Note 
that the expression of Cf. given in (5.20) for k € {0,1, . . . ,d - 1) follows from the second 
sum in (7.54). 



40 



MICHEL L. LAPIDUS, ERIN P. I. PEARSE, AND STEFFEN WINTER 



8. Concluding remarks and future directions 

8.1. Relation to previous results. We will now discuss in more detail the consistency of 
our tube formula with the tube formulas for fractal sprays and strings previously obtained; 
see also Remark 5.8. 

8.1.1. Comparison of the present pointwise results with the distributional results o/[LaPe3]. 
Recall that in [LaPe3], the tube formula obtained is only shown to hold distributionally, and 
only for fractals sprays with monophase generators. (For a discussion of how Theorem 4.1 
extends results of [LaPe3] to generators which may not be monophase (or even pluriphase), 
see Remark 2.3.) 

For monophase generators G, the tubular zeta function in Definition 3.5 simplifies 
to the zeta function appearing in [LaPe3, Definition 7.1], and consequently Corollary 4.2, 
the monophase case of Theorem 4.1, is precisely the pointwise analogue of [LaPe3, Theo- 
rem 7.4]. We leave this as an exercise to the reader, with the following hints: 

(i) Note that the constant Kd{G) has a different meaning in [LaPe3, Eq. (5.9)], namely 
Kii{G) - -Aj(G). In this paper, we have K^iG) = in the monophase case (cf. 
Remark 2.3) and Ad{G) is kept as /Irf(G) in the formulas. 

(ii) When one computes the residue of ^^(e, i) at 5 = li in the version of [LaPe3], a 
term appears which cancels the term Ad{G)i^jA.d) in (4.1). 

Note that for the earlier distributional results, the assumptions on the underlying fractal 
string £. are slightly weaker: the order of languidity is arbitrary, and fractal strings with 
D = d are permitted for fractal sprays in R''. The additional assumption D < d on the 
abscissa of convergence of ^ in Theorem 4.1 is necessary for the proof to hold and is 
similar to the assumption D < 1 in [LvF2, Theorem 8.7]. Note that one always has D < d 
for a fractal spray with finite total volume, as the latter is given by l^j.{d)AdiG). Although it 
is easy to construct a fractal spray with D = d, it follows from Proposition 5.1 that a self- 
similar tiling cannot satisfy D-d; indeed, this would violate the nontriviality condition. 

8.1.2. Comparison with the 1 -dimensional case. To see that the pointwise tube formula 
for fractal strings in [LvF2, Theorem 8.7] is a special case of our tube formula for fractal 
sprays in Theorem 4.1, let 7~ be a fractal spray in R, i.e., a geometric fractal string. Then 
the generator G is always a bounded open interval of length 2g (g being the inradius of G) 
and with a (monophase) Steiner-like representation V(G-e) - 2e, for < s < g, implying 
Ko{G) = 2 and K[(G) - 0. The fractal string X. = {^i, ^2, • ■ .) of the scaling ratios generating 
T corresponds to the fractal string X. = {A, ^2, • ■ ■) of the lengths tj :- 2g(j of the intervals 
used in [LvF2], whence (i{s) - (IgYi^jAs). Since we are in the monophase case, the 
tubular zeta function (j- simplifies to 

^ris,s)^s ^^(.)^_-_j.^^(.)-_, 

which is precisely the function appearing in [LvF2, Thm 8.7]. Moreover, the complex 
dimensions at which the residues are taken also coincide, except for the two integer di- 
mensions and 1 . However, the residue at 1 cancels for the same reasons as in hint (ii) 
above, and one can show that the residue at appears in the tube formula (4. 1) if and only 
if G W \ £)£(W) = W \ DiiW), just as in [LvF2, Thm 8.7]. Finally, we remark that, 
in the setting of geometric fractal strings, the hypotheses of Corollary 4.2 are exactly the 
same as in [LvF2, Thm 8.7]. 
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8.2. Origin of the terms in the tube formula. The proof of Corollary 5.9 (given in Sec- 
tion 7.2.2) explains the origin of each term in the exact tube formula (5.18). Indeed, in 
(7.53)-(7.54), the first and second sum express the contribution of the tail zeta function 
5riaii(^' ■ ) '^he scaling and integer dimensions of T, respectively, while the third sum ex- 
presses the contribution of the residues of the head zeta function ^7-i,ead(^' ' ) the integer 
dimensions. 

8.3. The monophase case. Note that if G is monophase, its coefficient functions ki.(G, e) 
are constant (and equal to KkiG)). Consequently, the functions fk in (7.1) vanish identically, 
and hence so does Crhr^A^^^ (7.48). As a result, one has = (.j-f^^f, which is the case 
treated in [LaPe3]. This is so, in particular, when d = I and G is a bounded interval 
(i.e., in the case of a fractal string). As a result, the contributions of the residues of the 
head tubular zeta function (J-^,^.^^ vanish identically and thus do not have to be taken into 
account. Note that in the monophase case, one must also have lime^o+ KdiG, s) = 0, and 
hence Kd{G,e) - for all < £ < see also the discussion of the monophase and 
pluriphase case in Remark 2.3. Consequently, this explains why the c/* term drops out of 
the fractal tube formula appearing in [LaPe3]. 

8.4. The general case. In Remark 4.4, it was observed that it is extremely useful to be 
able to drop the assumption that S (0) < 0, especially for investigating delicate questions 
concerning the Minkowski measurability of fractal sprays and self-similar tilings. Indeed, 
by analogy with [LvF2, Section 8.4], the proof of such a result requires screens lying 
arbitrarily close to the line Re s - D. If in the languid case of Theorem 4. 1 one drops the 
requirement that 5(0) < 0, then the tube formula in (4.1) still holds, provided the error 
term Kis) given in (4.2) (or equivalently, given by '??taii(e) in (7. 19)) is replaced by 

!R(e)=9?taii(e) + ^head(e), (8.1) 

with Khsadis) as in (7.18). However, while the estimate (9(e''"™P'^) as s — > 0^ remains true 
for ^taii(e)j it will not be satisfied in general for !??head(e)^ the sum of the residues of 
over the hidden integer dimensions. Hence, such a tube formula would be rather useless, as 
its error term may be of the same order as (or even dominate) its 'main term'. As a result, 
this generalization of Theorem 4. 1 would not be suitable for investigating the Minkowski 
measurability of fractal sprays or even of self-similar tilings. 

In fact, the assumption 5 (0) < should be seen as the price one has to pay for the gen- 
erality of the allowed Steiner-like representations. Stronger hypotheses on the generator G 
(or on the coefficients in the Steiner-like representation) will lead to better estimates of the 
error term and thus allow one to drop this assumption, as in the monophase case, and to 
extend the results on Minkowski measurability mentioned in Remark 4.4 (and discussed in 
detail in [LPW2]) beyond the monophase setting. We plan to address this issue in [LPWl]. 

8.5. Piecewise analytic Steiner-like representations. In Example 6. 1, there is a partition 
of the interval (0,^] into finitely many pieces, namely iQ,g] - {0,g/ V2j U (^g/ V2,gj, such 
that each coeflicient function kii(G, s) is analytic on the interior of each subinterval. That 
is, each K^iG, s) is continuous and given by an absolutely convergent power series in s (in 
the first subinterval) or ^ (in the remaining subintervals). In such a case, we say that G has 
a piecewise analytic Steiner-like representation. 

This condition appears to be satisfied by many natural examples of fractal sprays (and 
self-similar tilings in particular). Indeed, it may be the key assumption needed to be able 
to apply our tube formulas efficiently to a wide variety of examples. In some future work. 
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we plan to investigate this property further, especially with regard to associated Minkowski i 

measurability results; see [LPWl] and Section 8.4 just above. 2 

8.6. Fractal curvatures. As was mentioned at the end of Section 1.1, a key motivation 3 

for the present work is the search for a good notion of fractal curvature. (See [LvF2, Sec- 4 

tion 8.2 and Section 12.7] for a discussion in the 1-dimensional case.) In our context, 5 

this would entail obtaining a local tube formula (with or without error term) correspond- e 

ing to Theorem 4. 1 (and its corollaries). This would lead naturally to an interpretation of 7 

the coefficients of such a local tube formula in terms of "curvature measures" (or rather, a 

distributions) associated with each complex dimension (that is, with each scaling and in- 9 

teger dimension). We hope to explore such a possibility in future work and to establish in 10 

the process some useful connections with [Win] and some of the references (on geometric 11 

measure theory and differential geometry) discussed in Section 1.1; see [Fed, Sch2, HLW] 12 

and [Wey, BG, Gra], in particular. Furthermore, we expect that eventually the present work 13 

and its ramifications will be helpful in obtaining global and local tube formulas (and an ap- 14 

propriate notion of fractal curvature), for more general fractal objects than fractal sprays 15 

and self-similar tilings. 16 
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